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ABSTRAGT 

Eclipsing systems of massive stars allow one to explore the properties of their components in great 
detail. We perform a multi-wavelength, non-LTE analysis of the three components of the massive 
multiple system S Ori A, focusing on the fundamental stellar properties, stellar winds, and X-ray 
characteristics of the system. 

The primary’s distance-independent parameters turn out to be characteristic for its spectral type 
(09.5 H), but usage of the Hipparcos parallax yields surprisingly low values for the mass, radius, 
and luminosity. Gonsistent values follow only if S Ori lies at about twice the Hipparcos distance, 
in the vicinity of the cr-Orionis cluster. The primary and tertiary dominate the spectrum and leave 
the secondary only marginally detectable. We estimate the V-band magnitude difference between 
primary and secondary to be AV ~ 2U8. The inferred parameters suggest the secondary is an early 
B-type dwarf (rs B1 V), while the tertiary is an early B-type subgiant (~ BO IV). We find evidence 
for rapid turbulent velocities (~ 200kms“*) and wind inhomogeneities, partially optically thick, 
in the primary’s wind. The bulk of the X-ray emission likely emerges from the primary’s stellar 
wind (logLx/ABoi ~ —6.85), initiating close to the stellar surface at Rq ^ 1.1 i?*. Accounting for 
clumping, the mass-loss rate of the primary is found to be log M Ri —6.4 [Mq yr“*], which agrees with 
hydrodynamic predictions, and provides a consistent picture along the X-ray, UV, optical and radio 
spectral domains. 

Subject headings: stars: individual (i5 Ori A) — binaries: close — binaries: eclipsing — stars: early- 
type — stars: mass loss — X-rays: stars 


1. INTRODUCTION 

Massive stars (M > IOMq) bear a tremendous in¬ 
fluence on their host galaxies, owing to their strong 
ionizing radiation and powerful stella r win ds (e.g., 
iKudritzki fc Pulsl 120001 : iHamann et al.l l2006f ). Yet 
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our understanding of massive stars and their evolu¬ 
tion still leaves much to be desired: (1) Values of 
mass-loss rates derived in different studies may dis¬ 
agree with e ach other b y up to an order of mag¬ 
nitude le.g.. iPuls et al.l 1996t iFullerton et al.l l2006t 
IWaldron fc Cassinellil I201C : IBouret et al.ll2012H and of¬ 
ten do not agr ee with theoreti cally predicted mass- 
loss rates (e.g., I Vink et HI 1200(1(1 . (2) The extent of 

wind inhomogeneities, which greatly influence mass- 
loss rates inferred by means of spectral analyses, are 
still largely debated (|Shavivl l2000f : lOwocki et al.l 12004 
Oskinova et al.l[2007t iSundavist et al.l[2011l : lSurlan et al.l 

201, ‘111 . (3) The production mechanisms of X-ray ra¬ 
diation in massive stars have been a central sub¬ 
ject of s tudy in recent decades (|Feldmeier et al.1 [19971 : 
iPollockI l20f)7fl and are still far from being under¬ 
stood. (4) The effect of magne tic fields and stellar 


rotation on massive stars (e.g.. iFriend & MacGregoi 

1984 

; (Maheswaran & Gassinelli 120091: Oskinova et al. 

2011 

:lde Mink et al.ll201.3l:IPetit et al.ll20l4 Shenar et al. 

2014 

), e.g., through magnetic braking (|Weber & Davis 

1967 

) or chemical mixing (|M^,(;^d 9 ^| 19 ? 7 f'- sDH being 


investigated. (5) Lastly, stellar multiplicity seems to play 
a fundamental role in the context of ma ssive stars, sig¬ 
nifica ntly affecting their evolution (e.g., lEldridge et al.l 

[Ml). 


Several studies in the past years le.g.. IMason et al. 

2000 

:IMafz ADellanidl2010l: IGhini et al.ll2012l: ISana et al. 

2013 

: iSota et al.112014 lAldoretta et al.l |2014(1 give di- 


rect evidence that at least half of the massive stars 
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are found in multiple systems. Massive stars in close 
binary systems generally evolve differently from sin¬ 
gle stars. S uch systems may experience si gnificant 
tidal for ces (iZahn [197^ iPaiate et al.l 1201,111 , mass- 
transfer (lPols et_aL Il991f l. additional supernova kicks 
( Hurley et al.l l 20 ?) 2 1 . and mutual irradiation effects 
I HowarthI Il997lb Given the large binary fraction, an 
understanding of these processes is critical to properly 
model the evolution of massive stars. Fortunately, bi¬ 
nary systems have two main advantages over single stars: 
First, they offer us the opportunity to empirically deter¬ 
mine stellar masses. Second, eclipsing binary systems 
provide the unique opportunity to investigate different 
characteristics o f stellar winds by taking advantage of 
occultation le.g.. lAntokhinl [ 201111 . It is therefore insight¬ 
ful to analyze eclipsing multiple systems of massive stars 
in our Galactic neighborhood, using adequate, state-of- 
the-art modeling tools. 

The star ,5 Ori A (HD 36486, Mintaka, HIP 25 930, 
HR 1852) is a massive triple system (see artist’s illus¬ 
tration in Fig.jT]) comprised of the close eclipsing binary 
Aa ( primary Aal, secondary Aa 2 ) with a 5.732d pe- 

riod (jHarvey et al.lll98^lHarvin et al.ll2?if)^ iMaver et al.l 

[20i3), and the more distant tertiary Ab at an an¬ 

gular separation of ~ 0 .3" relative to the binary Aa 
with a period of ^346 yr (iHeintdl 19801 : iPerrvman &: ESAl 
119971: iTokovinin et al.l 1201411 . The tertiary Ab has been 
photometrically resolved from the binary Aa in dif¬ 
ferent surveys (iHorch et "ni 120011 : iMason et al.l 120091 : 
iMaiz ADellaiiidl2010tl and is found to contribute ~ 25% 
to the system’s flux in the visual band. These three com¬ 
ponents are not to be confused with the more distant and 
significantly fainter stars <5 Ori B and C at separation 
33" and 53", respectively; together, these five stars com¬ 
prise the multiple system S Ori, also known as Mintaka. 
With a visual magnitude oi V = 2™24 outside eclipse 
([Morel fc Magnen^ll978ll . 5 Ori A is one of the bright¬ 
est massive multiple systems in the night sky, making it 
easily visible to the naked eye (it is the westmost star in 
Orion’s belt). Acco rding to the new Hipparcos reduction 
([van Leeuwenll200^ . its parallax is 4.71±0.58mas, corre¬ 
sponding to a distance of d = 212 ±30 pcQ. On the other 
hand, the system resides in the Orion OBlb association, 
to which the cr-Orionis cluster belongs as well. The dis¬ 
tance to the cluster itself is estimated to be d ~ 380 pc, 
almost a factor 2 larger th an the Hipparcos distance (see 
iGaballero fc Solan^l2008L and references therein). 

This paper is a contribution to a series of papers within 
the framework of the S Ori collaboration. The other pa¬ 
pers include the analysis of high quality Chandra obser¬ 
vations to explore X-ray properties (Corcoran et al. in 
prep. Paper I) and variability (Nichols et al. in prep. Pa¬ 
per H) in the system, and a complete photometric and 
spectrometric variabilty analysis (Pablo et al. in prep. 
Paper HI). In this study, we focus on a non-LTE, mul¬ 
tiwavelength spectral analysis of the three components 
Aal, Aa2, and Ab, with the goal of obtaining reliable 
stellar and wind parameters. 

5 Ori A has been repeatedly studied previously. The 


This value is n ot corrected for the Lutz- Kelker effect 
IILutz & Kelkeil 119731 1. IMafz Apellaniz et al.l II2008I 1 account for 
this effect and revise the distance to d = 221 pc. However, the 
difference is negligible within measurement uncertainties. 



Fig. 1. — Artist’s impression of the triple system <5 Ori A, as 
viewed from Earth. The primary (Aal) and Secondary (Aa2) form 
the tight eclipsing binary of period 5.7 d. The primary shows ev¬ 
idence for a significant wind in all spectral domains. A third star 
(Ab) at an angular separation of 0.3” (~ 100 AU at a distance 
of d = 380 pc to Earth) orbits the system with a period of ?a346 yr 
IIT'okovinin et ani2014l) and contributes roughly 25% to the total 
visual flux. The sizes of all three components are drawn are to 
scale, as inferred in this collaborative study. Their colors reflect 
their relative temperatures. Note that the tertiary is the second 
brightest companion in the system. The distance between the bi¬ 
nary system Aa and the tertiary Ab is not to scale. Credits to 
Jessica Mayo 


system, which shows clear evidence for a significant stel¬ 
lar wind in the optical, UV, and X- ray domains, ha s been 
assigned t he spectral type 0 9.5 H (|WalborDlll972l) . later 
refined bv iSota et al.l (I2011D to 09.5 H Nwk, which most 
probably corresponds to the brigtest component: Aal. 
With less confidence, the secondar y Aa2 has been as - 
signed the spectral type BO., 5 HI (iHarvin et al.l l2(i0^ . 
This result is questioned by IMaver et al.l ([201^ . who 
argue that the secondary is too faint for this spectral 
type. The latter authors further suggest the spectral 
type 09 IV for the tertiary, and leave the secondary un¬ 
classified. 

Although the stellar parameters of an eclipsing stel¬ 
lar system can usually be sharply constrained, much 
controve rsy is found in t he lit erature in the case of 
5 Ori A. iKoch fc Hrivn^ (I1981D reported Mi = 23 Mq, 
i?i = 17i ?0 and M 2 = 9Mq, R 2 = 10 i?© for the pri¬ 
mary and secondary masses and radii, respectively, as 
well as a V-band nr a gnitud e difference of AVAaiAa 2 = 
1±4. IHarvin et al.l (| 2002 f ) later inferred significantly 
smaller masses for the primary and secondary. Mi = 
11.2 M 0 , M 2 = 5.6 Mq, and a smaller contribution of 
the secondary in the visual band: AVAa,iAa 2 = 2±5 . 
These results were challenged by IMaver et al.l (|2010D . 
who suggested that a confusion between the secondary 








































































A non-LTE analysis of i5 Ori 


3 


Aa 2 and tertiary Ab led to the low masses obtained 
by iHarvin et all (|2002ll . iMaver et alJ (120101) inferred 
i?i = 15.6 i ?0 and i ?2 = 4.8 i?©. They did not detect 
any contribution from the secondary Aa2 and concluded 
that AVAaiAa 2 ^ 3™5. Assuming Mi = 25 M©, they in¬ 
ferred M 2 = 9.9 Mq. In Paper III, the secondary’s radial 
velocity (RV) curve could not be constructed, and, for an 
adopted primary mass of Mi Ri 24: Mq, it was concluded 
that M2 = 8.5 Mq Ri = 15.1 i?©, R2 = 5.0 Rq. 

Similarly, reported valu es of the mass-loss rate of the 
primary are quite diverse: iLamers fc Leithereil (jl993D re¬ 
port log Ml = —5.97 [M© yr“^] based on radio observa¬ 
tion and log Ml = —5.92 [M© yr“^] based on Ha analy¬ 
sis. ILamers et al.l (jl999D later similarly obtain log Mi = 
—5.7 [M© yr“^] based on P Cygni profile analysis using 
the Sobolev plus exact integration (SEI) method. Ne¬ 
glecting the effects of wind inhomogeneities and poros¬ 
ity (jOskinova et al.l[2007HSurlan et al.ll2013[ ) on the ob¬ 
served X-ray spectrum, and adopting a generic 0-type 
model, iCohen et al.l 1)20141) derived a value of Mi = 
—7.2 [M© yr“^] from X-ray line profile fitting, more than 
an order of magnitude less than the previously inferred 
value. Indeed, no consensus on the mass-loss rate of the 
primary has been reached so far. 

5 Ori A has an observed X-ray luminosity of 
Lx ~ 1.4 • 10^^ergs“^ for d = 380pc (Paper I). ^ 
Ori A’s X-ray prop er ties were previou sl y ex plored 
by iMiller et al.l (120021). Leutenegger et al.l (( 2 OO 6 D . and 
iRaassen fc Pollock 2013 ). These studies generally iden¬ 
tify the X-ray formation process to be intrinsic to the pri¬ 
mary’s wind, a result which is further supported within 
our collaboration. An extensive review of past X-ray ob¬ 
servations and analyses are given by Papers I and II. 

In this study, we perform a consistent non-LTE photo¬ 
sphere and wind analysis of the three components of the 
triple system 6 Ori A in the optical, UV, and X-ray do¬ 
mains, at several orbital phases. We analyze the optical 
and UV spectra using the non-LTE Potsdam Wolf-Rayet 
(PoWR) code, which is applicable to any hot star. We 
further illustrate the importance of optically thin and 
optically thick clumps in the wind. We use the non-LTE 
models to simulate the effect of X-rays in the wind of the 
primary and derive onset radii of X-ray formation regions 
using ratios of forbidden and intercombination lines in 
the Chandra spectra. Finally, using the non-LTE model 
of the primary, we calculate synthetic X-ray line profiles 
and compare them to observed ones. Our results are fur¬ 
ther compared to studies of the radio emission from the 
system. This study thus encompasses the whole range 
from the X-ray domain to the radio domain. 

The structure of the paper is as follows: In Sect. [5] 
we describe the observational data used. Our modeling 
methods and assumptions are discussed in Sect. [S] In 
Sect.m we present and discuss our results. Sect.[S]focuses 
on the effect of wind inhomogeneities on the optical, UV 
and X-ray spectral domains, while in Sect. |6]we study 
the X-ray radiation of the star. Lastly, in Sect. [71 we 
summarize our results. 


2. OBSERVATIONAL DATA 

All spectra used in this analysis contain the contri¬ 
bution of the primary Aal, secondary Aa2, and ter¬ 
tiary Ab. In the following, all phases given are pho¬ 


tometric phases relative to primary minimum (occurring 
when the secondary occults the primary), calculated with 
Eo = HJD 2456277.79 d and P = 5.732436 d (Paper III, 
and references therein). 

Two of our optical spectra are of high resolution and 
high signal-to-noise, obtained with the NARVAL spec- 
tropolarimeter on 23-24 October 2008. These spectra 
were reduced with sta ndard techniques a nd downloaded 
from the PolarBas43 ([Petit et al.ll20l4l . The observa¬ 
tions were carried out with the Telescope Bernard Lyot. 
The spectra have a S/N of 500 — 800, and correspond to 
phases (j) = 0.84 and (f) = 0.02. 

Three additional optical spectra at phases 
(j) = 0.19, 0.38, and 0.54 were obtained on the nights of 
28, 29, and 30 Dec 2012 (contemporaneous with out 
Chandra and MOST observations) using the CAFE 
spectrograph at the 2.2 m Calar Alto telescope as 


part of the CAFE-BEANS project, a survey that is 
obtaining high-resolution (i? ^ 65 000) multi-epoch 
optical spectroscopy of all bright O stars in the northern 
hemisphere to compleme nt OWN, the equiv a lent south¬ 
ern h emisphere survey (jBarba et al.l 120101 : iSota et all 
120141) . On each date, ten consecutive 30s exposures 
were obtained and combined to yield spectra with S/N 
varying from ~ 200 in the blue part (^ 4500 A) to ~ 500 
in the red part (^ 6000A). The data were processed 
using a pipeline developed specifically for the project. 
The velocity stability was checked using ISM lines. 
M ore details regarding the reduction pipeline are given 
bv iNegueruela et all (|2014D . 

For the spectral range 1200 — 2000 A, we make use of 
archival lUE spectra at different orbital pha ses (see the 
observation log given bv IHarvin et al.ll2002D . The spec¬ 
tra have a signal-to-noise (S/N) ratio of ~ 10/pixel, and 
are averaged in bins of 0.1 A. The lUE spectra are flux- 
calibrated and are rectihed using the model continuum. 

In the spectral range 1000 — 1200 A, we make use of 
the Copernicus observation available in the Copernicus 
archive under the ID c025-001.u2. This observation con¬ 
sists of 48 co-added scans obtained between 21 and 24 
November 1972. 

For the spectral energy dis tribution (SEP), we mak e 
use of U, B, V, R, /, J, H, K ([Morel h Magnenatlll978l). 
WI^ (iCutri fc: et al.llSoT^ . and IRAS ( Moshir fc et al l 
I 1990 D photometry. We further use a low-resolution, flux 
calibrated optical spectrum kindly supplied to us by 
S. Fabrika and A. Valeev (priv. com.). The spectrum 
was obtained with the Russian BTA telescope using the 
SCORPIO focal reducer, on 31 Dec 2013 in the range 
3800 — 7200 A with a spectral resolution of 6.3 A, and 
corresponds to phase (j) = 0.42. A Hartmann mask was 
used to avoid saturation. 

The Chandra X-ray spectra used in Sect. [6] were taken 
with the BEG and MEG detectors for a total exposure 
time of 487.7 ks. The data are thoroughly described in 
Papers I and H. 


3. NON-LTE PHOTOSPHERE AND WIND MODELING 
3.1. The PoWR code 


http://polarbase.irap.omp.eu 



























































4 


PoWR is a non-LTE model atmosphere code especially 
suitable for hot stars with expanding atmospheres^- The 
code consists of two main parts. In the first part, referred 
to as the non-LTE iteration, the co-moving frame radia¬ 
tive transfer in spherical symmetry and the statistical 
balance equations are solved in an iterative scheme un¬ 
der the constraint of energy conservation, yielding the 
occupation numbers in the photosphere and wind. The 
second part, referred to as the formal integration, deliv¬ 
ers a synthetic spectrum in the observer’s frame. The 
pre- specified wind velo city field takes the form of a /3- 
law (jCastor et al.lll9'^ in the supersonic region. In the 
subsonic region, the velocity field is defined so that a hy¬ 
drostatic density stratification is approached. Line blan¬ 
keting by the iron lines is t reated in the superlevel ap- 
nroach (iGrafener et al.l[ 200 ^ . as originally introduced by 
lAndersonT i 1989f l. A closer description of t he assumptions 
and in etho ds used in the code is given by IGrafener et al.l 
(|2002ll and iHamann fc Grafened ( 200411 . 

In the non-LTE iteration, line profiles are Gaussians 
with a constant Doppler width UDop- In the formal in¬ 
tegration, the Doppler velocity is composed of depth- 
dependent thermal motion and microturbulence. The 
microturbulence ^(r) is interpolated between the pho- 
tospheric microturbulence C(rph) = ^ph and the wind 
microturbulence ^(rw) = with the radii Cph and 
Cw pre-specified. Thermal and pressure broadening are 
accounted for in the formal integration. Turbulence 
pressure is also accounted for in the non-LTE itera¬ 
tion. Optically thin wind inhomogeneities are accounted 
for in the non-LTE i teration using the so -called ‘micro- 
clumping’ approach (|Hamann et al.l [199511 . The density 
contrast between a clumped and a smooth wind with 
an identical mass-loss rate is de scribed by the depth- 
depen dent clumping factor D{r) (jHamann fc Koesterkd 
119981) . where the clumps are assumed to be optically 
thin. Optically thick clumps, or ‘macroclumps‘, are ac- 
counted for in the fo rmal integration as described by 
lOskinoya et al.l (|2007f) . where the clump separation pa¬ 
rameter Linac is to be specified (see Sect. [5]). 

Four fundamental input parameters define a model at¬ 
mosphere of an OB type star: T*, L, g*, M. T* is the 
effective temperature of a star with luminosity L and 
radius i?*, as defined by the Stefan-Boltzmann relation 
L = AttctRIT^. 5 * is related to the radius R^, and mass 
M* via g* = g{R*) = G M^R~‘^. M is the mass-loss rate 
of the star. The stellar radius ii* is defined at the con¬ 
tinuum Rosseland optical depth rftoss=20, which is the 
inner boundary of our model, and which was tested to 
be sufficiently large. The outer boundary of the model 
is set to 100 A*. Note that the stellar radius is generally 
not identical to the photospheric radius R 2/3 defined at 
tross = 2/3. However, usually i?* = i? 2 / 3 ; except in cases 
of extreme radiation pressures (e.g., supergiants, Wolf- 
Rayet stars). Nevertheless, one must bear in mind that 
the effective temperature referring to the photospheric 
radius, which we denote with T 2/3 to avoid ambiguity, 
may slightly differ from T*. Similarly, 52/3 = ^(-^ 2 / 3 ) is 
the gravity at tross = 2/3. Like most studies, we specify 
photospheric values when compiling our results in Ta- 

PoWR models of Wolf-Rayet stars can be downloaded at 
http: //www.astro.physik.uni-potsdam.de/PoWR.html 


ble[Tl but the cautious reader should be aware of this 
difference when comparing with other studies. 

Due to the strong radiative pressures in massive stars, 
one cannot measure the gravity g* directly from their 
spectra, but rather the effective gravity ^eff- 9* is ob¬ 
tained from 5 eff via '■= g* (l — Frad), where F^ad is a 
weighted average of the ratio of total outward radiative 
force to inward gravitational force over the hydrostatic 
domain. The outward radiative force is calculated con¬ 
sistently in the non-LTE iteration, and includes the con¬ 
tribution from line, continnum, and Thomson opacities 
(Sander et al. submitted). 

3.2. The analysis method 

The analysis of stellar spectra with non-LTE model 
atmosphere codes is an iterative, computationally ex¬ 
pensive process, which involves a multitude of parame¬ 
ters. Nevertheless, most parameters affect the spectrum 
uniquely. Generally, the gravity g^ is inferred from the 
wings of prominent hydrogen lines. The stellar tempera¬ 
ture T* is obtained from the line ratios of ions belonging 
to the same element. Wind parameters such as M, Voo-, 
and Lniac are derived from emission lines, mainly in the 
UV. The luminosity L and the colour excess E{B — V) 
are determined from the observed spectral energy dis¬ 
tribution (SED), and from flux-calibrated spectra. The 
abundances are determined from the overall strengths of 
lines belonging to the respective elements. Finally, pa¬ 
rameters describing the various velocity fields in the pho¬ 
tosphere and wind (rotation, turbulence) are constrained 
from profile shapes and strengths. The radius i?* and 
spectroscopic mass M, follow from L, T* and g*. 

To analyze a multiple system such as S Ori A, a model 
for each component star is required. As opposed to single 
stars, the luminosities of the components influence their 
contribution to the overall flux and thus affect the nor¬ 
malized spectrum. The light ratios of the different com¬ 
ponents therefore become entangled with the fundamen¬ 
tal stellar parameters. Fortunately, with fixed tempera¬ 
tures and gravities for the secondary and tertiary compo¬ 
nents, the observational constraints provide us with the 
light ratios (see Sect. 13.31) . 

Methods to disentangle a composite spectrum into its 
constituent spectral components by observing the system 
at different phases have been pr oposed and implemented 
during the past few years (iBagnuo lo fc Gi e)^ 119911: 
Simon fc Sturml[l994lHadravalll995tlGonzalez fc Levatol 
2006l : lTorres et al.ll201lD . In fact. an attemp t to disentan¬ 
gle the spe ctrum of <5 Or i A wa s pursued bv iHarvin et al.l 
(|2002f) and iMaver et al.l (120101) . However, even after per¬ 
forming the disentanglement, the two sets of authors 
come to significantly different results. A disentanglement 
of the He I A6878 line was performed in Paper HI, but, 
after accounting for the contribution of the tertiary as 
obtained here, no clear signal from the secondary was 
detected. Given the very low contribution of the sec¬ 
ondary Aa2, and having only poor phase coverage in the 
optical, we do not pursue a disentanglement of S Ori A. 

As a first step, the secondary and tertiary models are 
kept fixed. Motivated by the results of Paper HI, we ini¬ 
tially adopt T 2/3 ^ 25 kK and log 52/3 4.15 [gcm“^] for 

the secondary. The tert iary is initia l ly fixe d with the 
parameters suggested bv iMaver et al.l (|2010D . With the 
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light ratios at hand, this hxes the luminosities of Aa2 and 
Ab. Having hxed the secondary and tertiary models, we 
turn to the second step, which is an accurate analysis of 
the primary. 

To constrain T*, we use mostly Hei and Hell lines, 
such as He I AA4026, 4144, 4388, 4713, 4922, 5015, 6678 
and Hell A4200, 4542, 5412, 6683. The prominent Hell 
A4686 lin e is found to be a poor temperature indicator 
(see Sect. 14.2|) . The temperature is further verihed from 
lines of carbon, nitrogen and silicon. 5 * is primarily de¬ 
rived from the wings of prominent Balmer and He ii lines. 
Here we encounter the difficult problem of identifying 
the contribution of the tertiary to the hydrogen lines 
due to the pronounced wings of the Balmer lines. We 
therefore also made use of diagnostic lines such as Cm 
A5696, whose behavior heavily depends on the gravity 
(see Sect. 14.41) . 

M and Tmac follow from a simultaneous htting of Ha 
and UV P Cygni lines. The wind parameters are checked 
for consistency with previous analyses of radio observa¬ 
tions and with X-ray observations (see Sect. 16.311 . The 
terminal velocity Voo is determined from UV resonance 
lines. However, Voo can only be determined accurately 
after the wind microturbulence has been deduced. The 
determination of the projected rotational velocity v sin i 
and of the inner (photospheric) and outer (wind) micro- 
turbulent velocities ^ph and ^w, as well as of the macro- 
turbulent velocity Umac, is discussed in detail in Sect. 14.31 

The abundances are determined from the overall 
strengths of lines belonging to the respective elements. 
We include the elements H, He, C, N, O, Mg, Al, Si, P, 
S, and elements belonging to the iron group (e.g., Fe, Ni, 
Cr etc.). With the remaining stellar parameters hxed, 
this is a straight-forward procedure. 

The total bolometric luminosity L = Li + L 2 + L 3 is 
obtained by htting the synthetic hux with the observed 
SED. The individual component luminosities follow from 
the light ratios. The color ex cess Eb-v is derived using 
extinction laws published bv iCardelli et al.l (1198911 with 
Rv = 3.1, and can be very accurately determined from 
hux-calibrated lUE observations. The inferred value for 
Ay is consistent with those of Maiz Apellaniz & Sota (in 
prep. 2015), who obtain a value of 0.185±0.013 mag using 
Tvcho-| -Johnson-|-2MASS photomet ry and the extinction 
laws of iMaiz Apellaniz et al.l (120141) . 

After obtaining a satisfactory model for the primary, 
we continue to iterate on the secondary and tertiary mod¬ 
els by identifying small deviations between the compos¬ 
ite synthetic and observed spe ctra at different orbital 
phases, as described in Sect. 14.21 This allows us to adjust 
the temperatures and gravities of both models, as well as 
of the projected rotation velocities. After improving the 
secondary and tertiary models, we return to the primary 
and adjust the model parameters accordingly. We re¬ 
peat this process several times until a satisfactory ht is 
obtained in the UV and optical, taking hundreds of lines 
into account in this process. 

3.3. Initial assumptions 

Given the large number of parameters involved in this 
analysis, it is advisable to initially hx parameters which 
are constrained based on observations, previous studies, 
and theoretical predictions. 


We adopt (3=0.8 for the exponent of the (3-law for all 
components, which is both supported by observations 
as w-ell as theoreticall y predicted (e.g., iKudritzki et al.l 
119891: iPnls et aD I1996D . Varying this parameter in the 
range 0.7 — 1.5 did not signihcantly effect on the result¬ 
ing synthetic spectrum. 

A V-band magnitude difference of AVAaAb = 1“35 
between the binary system Aa (composed of Aal and 
Aa2) and the t ertiary Ab was measure d by the Hiy- 
parcos satellite (|Perrvman fc ESTG I1997D . iHorch et al.l 
(| 200 H) report AVAaAh = 1 U59, obtained from speckle 
photometry. iMason et al.l (j2009[ l hnd AV AaAh = 1U4 
by me ans of speckle interferometry. Like lMayer et al] 
(iMoli . we adopt AVAaAb = 1™4, corresponding to the 
flux ratio RAaAb '■= Fv{Aa)/Fy{Ab) = 3.63. As for 
the binary components Aal and Aa 2 , additional infor¬ 
mation regarding AVAaiAa 2 can be obtained from the 
yisual light curye of the system S Ori A. The secondary 
light curve minimum (primary star in front) has a depth 
AVii.niin ~ 0U055 with an error of about 0U005. Since 
the secondary eclipse is partial (Paper HI), the secondary 
minimum yields a lower limit to RAa 2 Aai- One obtains 
after some algebra 


n ^ (^II .min l)(3?AaAb + 1) 

nAa2Aal > -^-—j- 

-^AaAb -an.min “r ^ 


( 1 ) 


In our case, i?AaAb ~ 3.63 and i?ii,min ~ 1.05, and 
so Fy{Aa2) > 0.07Fy(Aal), which in turn implies 
AVAaiAo 2 ^ 2U9. It follows that the secondary Aa2 con¬ 
tributes at least 6.5% to the visual flux of the binary, and 
at least 5.1% to the total visual flux of the system. We 
thus initially assume AV 4 a 2 Aoi = 2 U 8 , and further con¬ 
strain this value in the spectral analysis (see Section4.1). 

In Sect. 14.11 we show that usage of the Hipparcos dis¬ 
tance d = 212 pc results in extremely peculiar param¬ 
eters for the primary given its spectral type, and that 
a much better agreement is obtained for the alterna¬ 
tive distance of the neighboring stellar cluster g-Orionis, 
d ~ 380 pc (see also discussion by iSimon-Diaz et al.l 
l2015fl . We therefore refer to both distances in the follow¬ 
ing, and discuss this issue more thoroughly in Sect. 14.II 
We note, however, that our results can be easily scaled 
with the distance, should it be revised in future studies 
(see Sect.S]). 

As we illustrate in Sect. [5l there are several indica¬ 
tions for wind inhomogeneities (clumps) in the wind 
of Aal. While clumping may already initiate at sub- 
photosp heric layers e.g.. due to sub-photospheric con¬ 
vection ()Cantiello et al.l [20091) . we avoid an attempt to 
treat dumpiness in the optically thick layers below the 
photosphere. For reasons which are discussed in detail 
in Sect.[5l we adopt a depth-dependent clumping fac¬ 
tor D{r) which is hxed to 1 (no clumping) in the do¬ 
main i?* < r < 1.1 A* and reaches a maximum value of 
34max = 10 at r > 10 i?*. 

Due to the faintness of the secondary and tertiary, we 
we can only give upper limits to their mass-loss rates, and 
set their terminal velocities to be 2 .6 times their escape 
velocity (iLamers &: Cassinellilll999ll . We further assume 
^{r) = 10 kms“^ for the secondary, a typical value for 
OB type stars of luminosi ty classes HI and lower (cf. 
iVillamariz &: Herrerol 120001 and references therein). We 
choose to adopt an identical density contrast iA(?’) for all 
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TABLE 1 


Inferred stellar parameters for the 

MULTIPLE STELLAR SYSTEM S ORI A 


Aal 

Aa2 

Ab 

T 2/3 [kK] 

29.5±0.5 

25.6±3000 

28.4±1500 

logSeff [cms- 2 ] 

3.0±0.15 

3.7 

3.2±0.3 

log 92/3 [cms-2] 

3.37±0.15 

3.9 

3.5±0.3 

Voo [kms“^ ] 

2000±100 

1200 

2000 

D 

10 

10 

10 

Lmac [-R©] 

^0.5 

- 

- 

E{B - V) [mag] 

0.065±0.01 

0.065±0.01 

0.065±0.01 

Av [mag] 

0.201±0.03 

0.201±0.03 

0.201±0.03 

V sin 2 [kms~^ ] 

130±10 

150±50 

220±20 

?ph [kms-l ] 

20±5 

10 

10±5 

[kms“l ] 

200±100 

10 

10 

Vninc [kms“l ] 

60±30 

0 

50±30 

H (mass fraction) 

0.70±0.05 

0.7 

0.7 

He (mass fraction) 

0.29±0.05 

0.29 

0.29 

C (mass fraction) 

2.4± 1 X 10-3 

2.4 X 10-3 

2.4 X 10-3 

N (mass fraction) 

4.0±2 X 10“^ 

7.0 X 10-4 

7.0 X 10-4 

0 (mass fraction) 

6.0±2 X 10-3 

6.0 X 10-3 

6.0 X 10-3 

Mg (mass fraction) 

6.4 X 10-4 

6.4 X 10-4 

6.4 X 10-4 

A1 (mass fraction) 

5.6 X 10-5 

5.6 X 10-5 

5.6 X 10-5 

Si (mass fraction) 

4±2 X 10-4 

6.6 X 10-4 

6.6 X 10-4 

P (mass fraction) 

5.8 X 10-® 

5.8 X 10-® 

5.8 X 10-® 

S (mass fraction) 

3.0 X 10-4 

3.0 X 10-4 

3.0 X 10-4 

Fe (mass fraction) 

1.3 X 10-3 

1.3 X 10-3 

1.3 X 10-3 

Adopted distance d [pc] 

212 380 

212 380 

212 380 

logL [Lq] 

4.77±0.05 5.28±0.05 

3.7±0.2 4.2±0.2 

4.3±0.15 4.8±0.15 

logM [M 0 /yr] 

-6.8±0.15 -6.4±0.15 

-8.5 -8.1 

< -7.0 < -6.6 

Afspec [Af©] 

7.513.5 2411° 

2.6 8.4 

7.0I4 22.5I14 

R2/3 [.R©] 

9.2±0.5 16.5±1 

3.6±1 6.5II5 

5.8±1 10.4±2 

My [mag] 

-4.47±0.13 -5.74±0.13 

-1.7±0.5 -3.0±0.5 

-3.2±0.4 -4.5±0.4 


three models. Finally, we do not attempt to constrain 
the abundances in t he secondary and ter tiary, but rather 
adopt solar values (|AsDlund et al.11200^ . This set of as¬ 
sumptions should bear very little effect on the derived 
fundamental parameters of the three components. 

4. RESULTS 

The stellar parameters inferred for the three compo¬ 
nents are given in Table[T] The model of the primary 
includes the effect of macroclumping and X-rays, which 
we thoroughly discuss in Sects. [5] and The upper part 
of the table displays parameters which, to first order, do 
not depend on the adopted distance. The lower part of 
the table denotes distance-dependent parameters. For 
the convenience of the reader, we give these parame¬ 
ters for both “candidate” distances, d = 212 pc and 
d = 380 pc. Luminosities scale as L oc d^, mass-loss 
rates as M oc radii as i?* oc d, and the mass as 

Mspec oc d^. The error margins given in Table [1] are 
discussed in Sect. 14.41 and are based on sensitivity of 
our analysis to these parameters. Values without errors 
which are not upper bounds indicate adopted values. 

The upper panel of Fig. [5] displays the synthetic SEDs 
for the models of the three components of the system. 
The total synthetic SED (red solid line) comprises the 
synthetic fluxes of the primary Aal (green dashed line), 
secondary Aa2 (gray dotted-dashed line), and tertiary 
Ab (pink dotted line). Recall that the tertiary is brighter 
than the secondary. The blue squares denote U, B, 
V, R, I, J, H, K, WISE, and IRAS photometry. We 
also plot two observed flux-calibrated spectra in the UV 


and optical (blue lines). The lower panels of Fig. [2] show 
the composite synthetic (red dotted line) and observed 
(blue line) normalized UV and optical spectra of the sys¬ 
tem. The optical spectrum and the UV spectrum in the 
range 1200 — 2000A correspond to phase 4> ~ 0.54. The 
synthetic normalized spectrum consists of three models 
calculated for Aal, Aa2, and Ab, shifted according to 
their RVs. The synthetic spectra are convolved with a 
Gaussian of FWHM = 0.1 A to account for instrumen¬ 
tal broadening, inferred from fitting of the interstellar 
Nai AA5890, 5896.3 lines. We do not plot the individual 
spectra of the three components for clarity. Note that 
the wavy pattern seen in the observed spectrum (e.g., in 
the domain 5900 — 6500 A) is an artifact which originates 
in the connection of the echelle orders, not related to the 
stellar system. 

It is evident that both the synthetic SED and normal¬ 
ized spectrum agree well with the observed spectrum. 
A good balance is obtained for all Hei lines and Hen 
lines, as well as for the metal lines. The inferred param¬ 
eters for microturbulence, rotation and macroturbulence 
yield consistent line strengths and profiles over the whole 
spectral domain. The pseudo continuum formed by the 
iron forest in the range 1300 — 1800 A, as well as most 
photospheric and wind features, are well reproduced. 

The few features which are not reproduced very well 
are the Balmer lines, and especially H(5. has a sig- 
nihcantly smaller observed EW (~ 1.85 A) compared to 
the synthetic spectrum (rs 2.35 A). In fact, the observed 
EW of H(5 is somewhat smaller than typical for similar 









7 


A non-LTE analysis of i5 Ori 



1100 1200 1300 1400 1500 1600 1700 1800 1900 




k/A 


Fig. 2.— Upper panel: Comparison between observations (blue squares and lines) and the total synthetic SED (red solid line), which 
consists of the primary (green dashed line), secondary (gray dash-dotted line), and tertiary (pink dotted line) models. Lower panels: 
Comparison between the composite synthetic (red dotted line) and observed (blue line) normalized spectra in the UV and optical. Both 
spectra roughly correspond to phase 0 = 0.54. For clarity, we refrain from showing the contributions of each component model to the 
normalized spectrum. The wavy pattern seen in the observed spectrum (e.g., in 6000 — 6200 A) is an artifact caused by connecting the 
echelle orders, not related to the stellar system. 










































spectral types ('e.g.. iCananzi et ^Il993ll . and the ques¬ 
tion arises as to the cause. A significantly smaller grav¬ 
ity for the primary does not agree with the wing shape 
of the other Balmer and Hell lines and hardly affects 
the EW. Reducing the hydrogen abundance of the pri¬ 
mary implies a larger helium abundance, which is not 
consistent with the helium lines. Larger M values lead 
to very strong emission in Ha, which is not observed. 
The photospheric microturbulence has only negligible ef¬ 
fect on the EW of hydrogen lines. We therefore conclude 
that the Balmer lines are dilluted by the light of one or 
both of the other components. However, the low EW of 
Hd could only be reproduced when assuming very pecu¬ 
liar parameters for the tertiary, e.g., a very weak gravity 
(logg 2/3 ^ 3.2[gcm“^]), or very large mass-loss rates 
(logM > —5.5 [M 0 /yr]). Such parameters are not only 
hard to justify physically, but also not consistent with 
the remaining spectral lines. The problem is also seen, 
albeit to a lesser extent, in the lines Hy and H/3. Future 
observations should shed light on this peculiarity. 

4.1. Which distance is the correct one? 

It is reassuring that those fundamental stellar param¬ 
eters of the primary which do not depend on the dis¬ 
tance match w ell with its sp e ctral type. Interpolating 
calibrations by iMartins et al.l (j2005l l to an 09.5 H class 
yields T 2/3 = 29.3 kK and log 32/3 = 3.35 [gcm“^], which 
agrees with our results within the error margins. The pri¬ 
mary’s nitrogen and silicon abundances are f ound to be 
slightl y subsolar. It is interesting to note that lSota et al.l 
(|2011ll recently added the suffix “Nwk” to the spectral 
type of Aal, implying relatively weak nitrogen lines, 
which we confirm independently. 

For d = 380 pc, the distance- dependent par a meter s 
agree well with calibrations by IMartins et al.l (j2005[l . 
Interpolating t o lum inosity class H, calibrations by 
IMartins et al.l (|2005ll imply logL/L© = 5.3, M = 
25 Mq, i? 2/3 = 18.2 Rq, and My = -5“73 for an 09.5 H 
star, which is consistent with our results. Furthermore, 
this distance implies radii and masses for the primary 
and secondary which agree well with the results obtained 
independently in Paper HI. However, very peculiar val¬ 
ues are obtained when using the Hipparcos distance of 
d = 212 pc (see lower part of Tablc[T|). In fact, all three 
components appear to be peculiar when adopting the 
Hipparcos distance. While a membership in a close bi¬ 
nary system offers some room for deviations from stan¬ 
dard values, this would not explain why the distant ter¬ 
tiary Ab should be peculiar too. The fact that distance- 
independent parameters are not unusual raises even more 
suspicion regarding the Hipparcos distance. 

A similar di s crepa ncy is observed other bright stars: 
iHummel et al.l (|2013[ 1 analyzed the binary system C Ori- 
onis A (y = 1.79) and inferred a distance of d = 294 ± 
21 pc based on an orbital analysis, and d = 387 ± 54 pc 
based on a photometric estimate. Both these distances 
are significantly larger than the corresponding Hipparcos 
distance of d = 225^*^27 PC ijvan LeeuwenI [2?)07h . which 
encouraged them to discard the Hipparcos distance in 
their study. Another bright star for which the Hippar¬ 
cos parallax implies very peculiar stellar parameters is 
the prototypical O supergiant ( Puppis {V = 2.25); Its 
distance was suggested to be at least twice as large as 


implied by its measured p arallax (jPauldrach et al.ll2012l : 
iHowarth fc Stevenil2014ll . 

Returning to 6 Ori A, our estimate of d ^ 380 pc 
would correspond to a parallax of tt ^ 2.6 mas, deviating 
by ~ 3.5(7 from th e newly reduced H ipparcos parallax 
of 4.71 ± 0.58 mas (jvan LeeuwenI [2f)(T7l l . The probabil¬ 
ity for such a deviation to arise randomly is thus ex¬ 
tremely small (< 0.1%). Ground-based parallaxes for 5 
Ori A generally suffer from much larger uncertainties. 
For example, the Yale catalog gives tt = 9.7 ± 6.7mas 
(jvan Altena et all 1199511 . It is interesting to note that 
the origi nal Hipparc os catalog gives a parallax of 3.56 ± 
0.83 mas (|ESAlll99'^ . In this case, the deviation could be 
plausibly explained as a random error, with 2.6 mas being 
^ la away. While stellar multiplicity has been suggested 
to cause systern atic errors in parallax measurements (e.g. 
lSzabado'^ll997[ l. the argument is unlikely to hold here 
given the different timescale of the involved orbits com¬ 
pared to that of the parallax measurement (which is to 
say, a year). However, large errors in the parallax mea¬ 
surement are also expected to occur in the case of very 
bright stars, which could lead to saturation, resulting 
in an inaccurate estimate of the barycenter of the point 
spread function. Brightness is indeed a property S Ori A 
shares with the objects mentioned above. It is beyond 
the scope of the paper to judge whether the new Hip¬ 
parcos reduction suffers from underestimated errors, but 
the fact that several bright stars show a similar pattern 
should encourage thorough studies on the matter. 

The distance to S Orionis is difficult to directly mea¬ 
sure with modern techniques. The system is probably 
too bright for ground-based (e.g. RECON^^) as well 
as space-based instruments (e.g. the Hubble Space Tele¬ 
scope’s Fine Guidance Sensor; the Gaia mission). Our 
best hope is to obtain the RV curve of the secondary, 
which requires optical spectroscopy with a S/N > 1000 
(Paper HI). Then, with a direct measure of the system 
parameters, follow-up long-baseline interferometry may 
provide the angular extent of the orbit, allowing for an 
or bital parallax measur ement similar to that of C Ori 
A (|Hummel et al.l[2013fl . At 380pc, the angular extent 
of the orbit will be on the order of 0.5 mas and thus 
difficult to measure, as the smallest an gular separation 
for a binary yet resolved is 1.225 mas (iRaghavan et al.l 
l2009f) . However, this was done in the K— band with the 
CHARA Array, so the resolution in the i?—band should 
be sufficient to resolve the binary. Moreover, along with 
a light-curve analysis, and possessing both RV curves, 
one could accurately derive the masses and radii of both 
components, which, upon comparison with a distance- 
dependent spectral analysis, would supply another con¬ 
straint on the distance. For now, we suggest that the 
Hipparcos distance may be underestimated, but leave the 
question open for future studies to resolve. 

4.2. Constraining the parameters of Aa2 and Ab 

The five columns of Fig.j^ depict six prominent 
Hell and Hei lines, from left to right: HeiiA4686, 
HeiA4922, Hen A5411, HeiA5876, and the two adja¬ 
cent lines HeiA6678 and HeiiA6683. Each row de¬ 
picts a different orbital phase, from top to bottom: 
4> = 0.02,0.19,0.38,0.54,0.84. This time, we explicitly 

17 


WWW. recons. org 












































9 


He 


A non-LTE analysis of i5 Ori 

He I Hell 


0 = 0.02 


(p = 0.19 


(p = 0.38 


.. 

•\ 
•\ 
•v 
• V 

- . 

if 

it 

it 

f 

1$ 

ft 

0 

. 

4680 4685 4690 

xik 

' ' ' 1 

' ' ' ' 1 

: \ / 

As- f 

• 

\*/ 

1 1 1 1 1 1 

i:/ : 

•r 

7 

4680 4685 4690 

' ' ' ' ' 

~i 1 1 1 1 r 

' T 

IS 

it 

it 




- - 

/j 

^' 1 ' ' ' 


xiK 


5405 5410 5415 

XiA 



He I 


He I + He 



-.'.' '' 





V 


7 V 

II 

7 

- V 

If 

// 

( 


/ f 

; 

~.1 ' \ 


' 1 ' ' r 


X/A 


6675 6685 

X/A 




4680 4685 4690 


4920 4925 




6675 6685 


0 = 0.54 





5405 5410 5415 



0 = 0.84 






6675 6685 


Fig. 3.— The 25 panels show the contribution of the primary (green dashed line), secondary (gray solid line), and tertiary (pink dashed- 
dotted line) models to the composite synthetic spectrum (red dotted line) for six prominent helium lines and five different phases, compared 
with observations (blue line). The lines depicted are, from left to right. He II A4686, He I A4922, He II A5411, He I A5876, and the two adjacent 
lines HelA6678 and HellA6683. The phases are, from top to bottom, (p = 0.02,0.19,0.38,0.54, and 0.84. At each phase, the component 
models are shifted with the velocities given in Table [2] 


show the relative contributions of the primary (green 
dashed line), secondary (gray solid line) and tertiary 
(pink dotted-dashed line) to the total synthetic spectrum 
(red dotted line), compared to the observations at each 
phase (blue line). In Table [21 we specify the RVs with 
which the three components are shifted at each phase. 
The RVs for the primary were inferr ed in this study and 
agree very well with the RV curve of iMaver et al.l (|201ClD 
and those of Paper III. 

4.2.1. The tertiary 

There is a wide and shallow spectral feature which does 
not originate in the primary and which is co nstant along 
all orbital phases (see Figs. [3] and [4]). Like iMaver et abl 


dMa), we identify this feature with the tertiary Ab, 
which is in fact the second brightest source in the system. 
The RV of the tertiary, which is practically constant over 
all phases, is also inferr ed independently, a nd agrees well 
with that suggested bv IMaver et al.l (I2010D . 

We find that the tertiary contributes significantly more 
to Hi lines than to Hen lines. Together with the vi¬ 
sual flux ratio implied from observations, this leads to a 
tertiary temperature of T* ~ 29.0 kK. The gravity and 
mass-loss of the tertiary were constrained based on the 
Balmer lines. As already discussed, it is very hard to 
identify the explicit contribution of the tertiary to these 
lines, and so the gravity and mass-loss of the tertiary are 
only roughly constrained. The parameters derived for 
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TABLE 2 

Radial velocities derived/adopted for each phase in 

KM 


Phase 

Aal^ 

Aa2A 

Ab^ 

0.02 

0 

10 

25 

0.19 

-92 

230 

25 

0.38 

-22 

120 

25 

0.54 

83 

-100 

25 

0.84 

105 

-210 

25 


'‘Values for primary and tertiary derived in the analysis 
'‘Values for secondary adopted from IHarvin et aTI KOO^ l. with the 
exception of phase 0.84, which is derived here. 

the tertiary (cf. table [T]) are consistent with it being a 
BO ly type star (IHabete fc Heintz^ll981l : ISchmidt-Kalerl 

fT^ 


are not observed), we are forced to change other stel¬ 
lar parameters, e.g., T, and usint. The secondary’s 
projected rotation velocity (tisint « 150kms“^) agrees 
with the feature shown in the right panel of Fig. IH and 
the secondary temperature (T Ri 26 kK) is consistent 
with that obtained from the light curve analysis of the 
system (P aper III). 

In Sect. 13.31 we argued that the light curve of S Ori Aa 
implies that the secondary contributes at least 5.4% to 
total visual flux of the system. In this section, we ar¬ 
gue that the secondary can contribute no more than this 
amount. In other words, our lower bound for the rel¬ 
ative flux contribution becomes also our upper bound, 
and therefore AVAaiAa 2 ~ 2™8. Together with the 
temperature of the secondary, this enables us to in¬ 
fer an approximate value for the luminosity of the sec- 
ondary. Its parameters suggest it is a Ri B1 V type star 
(|Habets &: Heintz3ll981t ISchmidt-KaledflOS^ . 


4.2.2. The secondary 




Fig. 4.— The two panels depict two observed Hei lines (blue 
lines) at phase (f) = 0.84 in which the secondary may be detectable. 
The different curves correspond to the primary (green dashed line), 
secondary (gray solid line), tertiary (pink dashed-dotted line), and 
total (red dotted line) synthetic spectra. 


It is very hard, or perhaps impossible, to recognize 
any contribution from Aa2 to the spectrum. One ex¬ 
ception might be the He i lines AA4026,4144 at phase 
(j) = 0.84, which we show in the left and right panels of 
Fig. m respectively (colors and lines are as in Fig. [I. At 
this phase, the primary’s RV approaches its maximum 
of ^ II0kms“^, and so the secondary is expected to 
be more easily observed. The He i A4026 line, for exam¬ 
ple, seems to have an extended wing towards blueshifted 
wavelengths. The very weak Hei A4I44 line is one of the 
few spectral lines which possibly portrays a partially iso¬ 
lated feature originating in the secondary. We therefore 
infer a RV for the secondary in this phase. For all other 

J ihases, we adopt the secondary’s RVs from IHarvin et alJ 
200 ll . since the secondary’ s lines canno t be is olated in 
the spectrum. Howe ver, as iMaver et alJ (|20I0ll pointed 
out, it is likely that iHarvin et al.l 1 20^ confused the 
secondary with the tertiary, so that the secondary RVs 
reported by IHarvin et al.l (1200211 are questionable. 

Even without directly detecting the secondary, we can 
still constrain its stellar parameters. Since the light curve 
provides a lower limit for its visual flux (see Sect. 13.31) . the 
luminosity of Aa2 cannot be arbitrarily small. Instead, to 
avoid too-strong line features from the secondary (which 


4.3. Bulk and turbulent motions 
4.3.1. Rotation and macroturbulence 

Rotation is usually the dominant broadening mecha¬ 
nism of photospheric metal lines in OB type spectra, as 
is the case here. We infer the projected rotational ve¬ 
locity V sin i by convolving the synthetic spectrum with 
a rotation profile and comparing the line shapes to the 
observed spectrum. While flux convolution is a fair ap¬ 
proximation for photospheric spectra, it does not account 
for the effect of limb darkening and for the extended for¬ 
mation regions of some lines. Accounting for rotation in 
the formal integration is essential for a consistent inclu¬ 
sion of these effects, but is computationally expensive. 
We therefore use the convolution method to infer vsini, 
and, only after obtaining the best-fit value of Lisini, do 
we account for rotation directly in the form al integration 
(cf. iHillier et al.ll2012l : iShenar et al.ll201-il , assuming co¬ 
rotation up of the photosphere (tross =2/3) and angu¬ 
lar momentum conservation in the wind. Compared to 
simple flux convolution, the detailed treatment of rota¬ 
tion gene rally yields de eper lines with less elliptic pro¬ 
files ('e.g.. IUnsoldllI95^ and may significantly affect the 
inferred abundances, projected rotation velocities, and 
even fundamental stellar parameters. 

The effect of rotation on the spectrum is coupled 
to the effects of macroturbulence Vmar., which we 
model using a Radia l -Tang ential profile (|Gravl II975I : 
iSimon-Dfaz fc Herrerol 120071) . Macroturbulence does 
not enter the radiative transfer per definition. There¬ 
fore, like solid body rotation, macroturbulence con¬ 
serves the EWs of the lines. Typical val ues of Vr^ar ^ 
50kms~^are reported for OB type stars (iLefeyer et al.l 
l2007t iMarkova h Pulsl[20081 : iBouret et al.ll20I2h . While 
the origin of macroturbulence is not certain, it has been 
suggested to be a manifestation of collective pulsational 
broadening ()e.g.. lAerts et al.ll2009H . usinf and Umac are 
inferred simultaneously from the profile shapes of helium 
lines and metal lines. Not accounting for macroturbu¬ 
lence generally results in qualitatively different profiles 
from what are o bserved in photospheric lines (see e.g., 
e xample given b y iPulsII2008D . 

IHarvin et al.l (l2002f) inferred usini = 157 km s for 
the primary. While this value agrees well with the wings 
of prominent helium and metal lines, it leads to too-broad 





























































A non-LTE analysis of S Ori 


11 


Doppler cores. A significantly better fit is obtained with 
i^mac = 60kms“^ and usini = 130knis“^. For an incli¬ 
nation oil = 76° (Paper III) and a radius of i? = 16.5 Rq 
(table [T]), this could imply that the rotation period is ap¬ 
proximately synchronized with the orbital period of 5.7 d. 

For the tertiary, we find an optimal fit for v sin i = 
220kms~^ and = 5 0kms ~^, thus confirming the 
findings of iHarvin et al.l (|2002ll and iMaver et al.l (|201Clf) 
that the tertiary is a rapid rotator. For the secondary, 
motivated by the arguments discussed in Sect. 14.21 we 
estimate fsini ~ 150kms“^ . We adopt r^mac = 0 for 
the secondary, lacking any spectral lines from which it 
can be inferred. Any other typical values would bear no 
effect on our results. 

4.3.2. Microtubulence and terminal velocity 

In contrast to macroturbulence, microscopic turbulent 
motion enters directly into the process of radiative trans¬ 
fer and generally affects the EWs of spectral lines. In 
the non-LTE iteration, we do not specify the microtur¬ 
bulence explicitly, but rather the total Doppler width 
WDop of the opacity and emissivity profiles. tiDop thus 
determines the resolution of the frequency grid in the 
non-LTE iteration. This parameter generally has a neg¬ 
ligible effect on the obtained population numbers, with 
the extreme ex ception of the Hell A4686 line. As was al¬ 
ready noted by lEvans et al.l (j20n4f ) , this line reacts very 
strongly to changes in woop- We illustrate this in Fig. [Si 
where we show a segment of the optical spectrum con¬ 
taining the He ii A4686 line. In the figure, we depict three 
models corresponding to the primary calculated with pa¬ 
rameters identical to those given in Table [H but with 
I’Dop = 20, 40, and 60kms“^, respectively. It is evident 
that the He ii 4686 A line reacts remarkably strongly to 
WDop- The exact origins of this effect are still under inves¬ 
tigation (Shenar et al. in prep. 2015), but likely involve 
a feedback effect in the highly non-linear iterative solu¬ 
tion of the radiative transfer problem. The remaining 
Hell lines show a much weaker reaction in the opposite 
direction. Other lines hardly respond to changes of this 
parameter. This example shows that, overall, the choice 
of the parameter WDop is not critical for the fit, and that 
Hen A4686 is a poor temperature indicator. Based on 
this, WDop = 30km/s is used in the analysis, consistent 
with the inferred microturbulence (see below). 

In the formal integration, apart from including natural 
and pressure broadening, the Doppler width is separated 
into a thermal component Uth, which follows the temper¬ 
ature stratification in the model, and a depth-dependent 
microturbulence component ^(r), which is assumed to be 
identical for all ions. As described in Sect. 13.11 ^(r) is 
interpolated between the photospheric ^ph and wind 
turbulent velocities between two pre-specified radii, all 
of which are free parameters. 

Values of photospheric microturbulence reported for 
O giants range between ~ 10kms~^ and 30kms~^ (e.g., 
Gies &: LambertI [19921 iSmartt et al.l 1199^ iBouret et alJ 
2012t) . Since the photospheric microturbulence .^ph is 
rarely found to be larger than 30kms“^, its effect on 
the profile width is negligible compared to the effect of 
rotation. However, ^ph can have a very strong effect 
on the EW of spectral lines. The abundances are thus 
coupled to ^ph, and wrong turbulence values can eas- 



Fig. 5. — This figure illustrates the extreme sensitivity of the 
Hen A4686 line to the parameter UDop- Here we depict on^ the 
primary model, calculated with the parameters in Table ^ but 
with UDop = 20kms“^ (black solid line), 40kms“^(red dashed 
line), and 60kms“^ (green dotted line). The formal integration is 
performed with the same turbulent velocity as given in Table ^ 
Notice that most lines hardly react to this parameter. 

ily lead to a wrong e s timation of the abundances (e.g., 
IMcErlean et al.l 119981 IVillamariz fc Herrerol I2000f) . To 
disentangle the abundances and turbulence from, e.g., 
the temperature, we take advantage of the fact that 
different lines respond individually to changes in abun¬ 
dances, turbulence, and temperature, depending on their 
formation process. Fig. |6| shows an example. The left, 
middle, and right panels depict the Hei A5876 line as 
observed at phase cj) = 0.84 (blue line). In each panel, 
we show three different composite synthetic spectra (i.e. 
composing all three components), which were calculated 
with parameters identical to those given in Table [I] ex¬ 
cept for one stellar parameter of the primary. In the left 
panel, T* is set to 29.5, 30, and 30.5 kK. In the middle 
panel, the helium mass fraction is set to X^e of 0.2, 0.3, 
and 0.4. Lastly, in the right-most panel, we set ^ph to 
15kms“^, 20kms“^, and 25kms“^, respectively. The 
three composite spectra at the left and middle panels can 
hardly be distinguished from each other, portraying the 
insensitivity of the He i A5876 line to temperature and 
helium abundance. In the relevant parameter domain, it 
is mainly ^ph which influences the strength of the Hei 
A5878 line. By considering hundreds of lines at all avail¬ 
able orbital phases, we find that ^ph = 20 km s“^ provides 
the best results for the primary’s model. Similarly, we 
find that a microturbulence of 10kms“^ for the tertiary 
yields the best global fit. 

We find evidence for a rapid increase of the turbu¬ 
lent velocity in the primary right beyond the sonic point. 
The left panel of Fig. [7] shows the CiV AA 1548,1551 res¬ 
onance doublet, as observed in the lUE spectrum taken 
at phase cj) = 0.73. The observation shows a wide absorp¬ 
tion trough which extends to red-shifted wavelengths, 
and our task is to reproduce this feature. The black 
dashed line depicts the composite synthetic spectrum 
calculated with the parameters in Table[I] but without 
an increased wind turbulence in the primary model, i.e. 

= ^ph = 20kms“^ . The absorption trough is not 
reproduced. Increasing the terminal velocity only af¬ 
fects the blue edge of the line. Varying (3 in the domain 
0.7 — 1.5 does not lead to any notable changes in the 
spectrum. The only mechanism that is found to repro¬ 
duce the redshifted absorption trough is a rapid increase 
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Fig. 6.— Sensitivity of the Hei A5876 line to temperature, he¬ 
lium abundance, and microturbulence. The observed spectrum 
(blue line) at phase (p = 0.84 is plotted along with three composite 
synthetic spectra calculated with the parameters given in Tabled 
but with different temperature (left panel: T = 29.5, 30, 30.5 kK), 
helium abundance (middle panel: Xjje = 0.2, 0.3, 0.4), and photo- 
spheric microturbulence (right panel: = 15, 20, 25 km s“^ ) for 

the primary model (black dotted-dashed, red dashed, and green 
dotted lines, respectively). 



Fig. 7. — Left panel: observed Civ resonance doublet at phase 
(p = 0.73 (blue line) is compared to two synthetic composite 
spectra, with increased wind turbulence of = 200 km s~^ in 
the primary model (red solid line), and without it (black dashed 
line). Large microturbulent velocities beyond the sonic point are 
the only mechanism that can reproduce the redshifted absorp¬ 
tion trough. Note that a further improvement of the profile is 
achieved by accounting for macroclumps, which are not included 
here. Right panel: Comparison of three lUE observations at phases 
p = 0.18,0.73,0.98 (black, red, and green lines, respectively), il¬ 
lustrating that the redshifted absorption trough is observed at all 
phases. 

of the microturbulence beyond the sonic point. The red 
solid line was calculated like the dashed black line, but 
with a wind turbulent velocity of = 200kms“^. ^{r) 
is assumed to grow from ^ph at rin < 1.1 R* to at 
J’out > 2 i?*. At the same time, the blue absorption edge 
is shifted by ^ 200 kms“^, thus influencing the value de¬ 
duced for Vac- The right panel of Fig. [7] shows the same 
CIV resonance doublet, as observed in three lUE spectra 
taken at phases ij) = 0.18,0.72,0.98 (black, red, and green 
lines, respectively). The figure illustrates the relatively 
small variability of this line, showing that our results do 
not depend on phase, and rejecting the contamination by 
another component as an explanation for the extended 
red-shifted absorption. The absorption trough is not re¬ 
produced for Tout significantly larger than 2 i?*, which 
is understandable given the need for red-shifted absorp¬ 
tion. It is interesting to note that it is microturbulence, 
and not macroturbulence, which is needed to reproduce 
this feature. A further improvement of the line profile 
fit is obtained by accounting for optically thick clumps 
(macroclumps) in the wind, as we will discuss in Sect. O 
We do not include macrocluinping at this stage in order 
to single out the effect of on the line prohle. 

Having inferred the turbulent velocity and after ac¬ 


counting for clumping in the wind of the primary, it is 
straight-forward to derive the terminal velocity Voo from 
resonance P Cygni lines. All prominent lines in the UV 
imply the same value for Voo (2000 ± 100kms“^). 

4.4. Uncertainties 

Since the calculation of a PoWR model atmosphere 
is computationally expensive, a statistical approach for 
error determination is virtually impossible. The errors 
given in Table[T] for each physical quantity are obtained 
by fixing all parameters but one and varying this param¬ 
eter, estimating upper and lower limits that signihcantly 
change the quality of the fit in many prominent lines 
relative to the available S/N ratio. Errors for parame¬ 
ters which are implied from fundamental parameters via 
analytical relations, e.g., the spectroscopic mass, are cal¬ 
culated by means of error propagation. In the case of 
multiple systems, all models but one test model are kept 
fixed, and only one parameter of the test model is varied. 



Xlk x/k 

Fig. 8.— Left panel: Sensitivity of Cm A5696 to gravity. The 
observed spectrum at phase (p = 0.02 (blue line) is compared to 
three composite synthetic spectra calculated with the parameters 
given in Table[T] but with logg* of the primary set to 3.9, 3.55, 
and 3.2 (black dotted-dashed, red dashed, and green dotted lines, 
respectively). Right panel: Sensitivity of He I A6678 to tempera¬ 
ture. We set the temperature of the primary to 29.5, 30.0, and 
30.5 kK (black dotted-dashed, red dashed, and green dotted lines, 
respectively). 

As mentioned in Sect. 21 it is very hard to constrain the 
gravity of the primary due to contamination by the other 
components. However, we took advantage of the fact 
that specific lines do drastically change their strength as 
a function of gravity. An example for such a diagnos¬ 
tic line, C III A5696, is shown in the left panel of Fig.[8l 
as observed at phase </) = 0.02 (blue line). Three com¬ 
posite spectra (i.e. including all components) are plot¬ 
ted, where only the logi/* of the primary is changed to 
3.9 (black dotted-dashed line), 3.55 (red dashed line), 
and 3.2 [gcm“^] (green dotted line). The remaining stel¬ 
lar parameters are kept fixed to the values given in Ta- 
blelll This line only starts to portray emission when 
logg* ^ 3.5 [gcm“^] . This line also serves as a good ex¬ 
ample to why fitting would not always suggest the best 
fitting model. The contribution of such a small line to 
the reduced is negligible, unlike its diagnostic power. 
The right panel of Fig.[5]depicts the sensitivity of the He i 
line to the stellar temperature. The temperature of the 
primary is changed to 29.5 (black dotted-dashed line), 
30 (red dashed line), and 30.5kK (green dotted line). 
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Most He I and He ii lines react strongly to changes in the 
temperature and thus enable us to sharply constrain it. 

5. INHOMOGENEITIES IN THE PRIMARY’S WIND 

Evidence for wind inhoinogeneities (clumping) in the 
winds of h ot star s are frequently reported. iHillied 
(|1984f) and iHillied (119911) illustrated the effect of op¬ 
tically thin clumps on the electron scattering wings 
of Wolf-Rayet emission lines. More compelling di¬ 
rect evidence for clumping in the form of stochas¬ 
tic variability on short timescales was observ ed in 
both Wol f-Rayet (e.g.. iLepine fc MoffaH I199E ) and 
OB stars (lEversberg et al.l Il998t iMarkova et al. 120051 : 
iPrinia fc Massall2010t ). Clumn sizes likely follow a con¬ 
tinuous distribution (e.g., power-law) which is intimately 
conn ected with th e turbulence prevailing in the wind 
(e.g.. iMoffati 11993) . However, since consistent non-LTE 
modeling of inhomogeneous stellar winds in 3D is still 
beyond reach, the treatment of clumping is limited to 
approximate approaches. 

5.1. Microclumping 

A systematic treat ment o f optic ally thin clumps 
was introduced by IHillied (jl984 ) and later by 
iHamann fc Koesterkd (| 19981 ) using the so-called micro- 
clumping approach, where the population numbers are 
calculated in clumps which are a factor of D denser than 
the equivalent smooth wind. In this approach, processes 
sensitive to p, such as resonance and electron scatter¬ 
ing, are not sensitive to clumping, and imply M di¬ 
rectly. However, processes which are sensitive to such 
as recombination and free-free emission, are affected by 
clumping, and in fact only imply the value of M\/13. 
This enables consistent mass-loss estimations from both 
types of processes, and offers a method to quantify the 
degree of inhomogeneity in the wind in the optically thin 
limit. 

To investigate wind inhomogeneities in the primary 
Aal, we first adopt a smooth model. The four panels of 
Fig. [9] depict four “wind” lines: three UV resonance dou¬ 
blets belonging to Siiv, Pv, and CiV (oc p), and Ha, as 
the only recombination line potentially showing signs for 
emission (oc p^). In each of the four panels, four com¬ 
posite synthetic spectra (i.e. containing all three com¬ 
ponents) are plotted along the observation (blue line). 
The four models differ only in the mass-loss rate of Aal: 
logM =-5.6 (black dashed line), -6.0 (red solid line), -7.1 
(green dotted line), and - 8 . 6 [M 0 yr“^] (purple dotted- 
dashed line). The remaining stellar parameters are iden¬ 
tical to the ones given in Table [T] (for d = 380 pc), but 
with D = 1 for the primary. The observations roughly 
correspond to phase (j) ~ 0.8, except in the case of the P v 
doublet, which is obtained from the co-added Copernicus 
data (see Sect. E]). 

While a best fit for Ha (lower right panel in Fig.jS]) 
is obtained with logM Ri — 6.0 [MQyr“^], we note that 
Ha portrays noticeable variability with time, which in¬ 
creases the uncertainty of the inferred mass-loss rate. 
Furthermore, our estimation of the mass-loss rate from 
Ha could be inaccurate due to underestimated contribu¬ 
tion to the emission by either the tertiary or by wind- 
wind collision effect. Fortunately, the Ha mass-loss rate 
is supported by radio observations. Assuming a smooth 


wind, iLamers fc Leithered (I1993D infer a mass-loss rate of 
logM Ri —5.95 [MQyr“^] based on the radio flux of the 
system for an adopted distance oi d = 500 pc. This value 
is revised to log M r; —6.1 [Mq yr“^] for d = 380 pc, only 
slightly smaller than the value derived for Ha in this 
study. Since both these processes scale as p^, and since 
the wind is assumed to be smooth in both cases, we there¬ 
fore conclude that log r; —6.0 [Mq yr“^] for the 

primary. 
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Fig. 9.— Observed “wind” lines (blue line): The Siiv 
AA1394,1408 (upper left), Civ AA1548,1551 (lower left), and Pv 
AA1118,1128 (upper right) resonance doublets, and Ha (lower 
right), roughly at phase (p ~ 0.8, except for the Copernicus 
data (P v) which are co-added. Each panel depicts four compos¬ 
ite spectra with the same parameters as in table ^ but without 
clumping {D = 1), and with different mass-loss rates: logM = 
—5.6, —6.0, —7.1, and -8.6 [Mq yr“^] (black dashed, red solid, green 
dotted , and purple dashed-dotted lines, respectively). The Ha and 
different P Cygni lines clearly imply different mass-loss rates, and 
cannot be fitted simultaneously. 


We now turn to the UV resonance lines. The Si iv line, 
shown in the upper left panel of Fig.[9l is clearly not 
saturated in the observation. In the model, it remains 
saturated for logM > — 6.8 [MQyr“^]. A best fit is 
obtained with logM —7.5 [Mq yr“^]. Together with 
the condition discussed in the previous paragraph, this 
implies D > 40; a best fit is obtained for D ~ 1000. The 
P V resonance line (upper right panel) implies a similar 
mass-loss rate, and hence a similar clumping contrast. 
The Civ line (bottom left panel), which does not look 
saturated in the observatioiF^ , requires a much lower 
mass-loss rate, of the order of logM — 8.0 [MQyr“^], 

We note that the lUE observation may suffer from calibration 
problems, which sheds doubt on whether the Civ resonance line is 
unsaturated. However, the observed shape of its absorption trough 
suggests it is indeed unsaturated. This is evident in each of the 60 
available lUE spectra, some even showing this more extremely. 
While we speculate that the line is desaturated here, our results 
do not depend on this strongly. 
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and as a consequence implies D k. 10^. 

But are such large density contrasts physically sound? 
ID and 2D time-dependent hydrodynamic simulations 
of line driven winds suggest typical values of I? = 
4 — 10, with D = 20 — 100 occurring in the most ex¬ 
treme cases (lOwocki et all 19881 iFeldmeier et al. Il997t 


iRunacres fc Owockil 1200^ ISundovist fc Owockil I2013I1 . 

Evolutionary considerations support D ~ 2 — 3 in or¬ 
der t o obtain suffi cient mass-loss from OB type stars 
(e.g., iHirschil [20081 1 ■ Values of 4 — 10 are typically re¬ 
ported for Wolf-Rayet stars based on electron scatter¬ 
ing wings of strong emission lines le.g.. lTodt et alJ[2013t 
iHainich et al.l[20T41 , while larger density contrasts of the 
order of 20 or mo re are suggested for OB stars (e.g., 
iBouret et al.ll201§l to reconcile the strong discrepancy 


between M values derived from P Cygni resonance lines 
and recombination lines. Analyses of the P v resonance 
doub let simultaneously wit h optical recombination lines 
(e.g.. iFullerton et'aril2006ll may imply clumping factors 
as large as 100. However, we no te that the latter au - 
thors neglect the effect of porosity (jOskinova et ahll^OOTH . 
and that the resulting low mass-loss rates are inconsis¬ 
tent wTHi__golnjizatmn_studles based on electron scatter¬ 
ing (|St-Louis fc Moffa'tll2008[l , which depends linearly on 
density. It seems therefore that the implied values for D 
exceed plausible limits, and that no single value for D 
can satisfy all emission features simultaneously. 

One possibility is that the abundances, or any other 
fundamental parameters derived for Aal which may af¬ 
fect the formation of the UV resonance lines, are signif¬ 
icantly inaccurate. This is unlikely, because the photo- 
spheric features, as well as the general shape of the UV 
iron forest, are well reproduced and do not leave room 
for variations. Another possibility is that one of the com¬ 
ponents, perhaps through interaction with the primary’s 
wind, is contaminating the observation. However, there 
are no clear indications for a periodic variation of the 
UV resonance lines with phase, as could be expected 
in such a case. We therefore suggest that not only is 
the wind of Aal clumped, but also that the clumps are 
op tically thick in t he st rong UV lines. As illustrated 
bv lOskinova et all (j2007f l. the optically thickness of the 
clumps leads to an effective reduction of the opacity of 
strong UV resonance lines, and could thus enable to ob¬ 
tain results which are consistent over the X-ray, UV, op¬ 
tical, IR and radio regimes. A further indication for op¬ 
tically thick clumps comes from the variable EW ratio of 
the of the two Si iv components. A ratio of 1 is expected 
to occur in the extreme case where the porous wind con¬ 
sists of optically thick clumps and “holes” between them, 
while a ratio of ~ 2 (corresponding to the ratio of oscilla¬ 
tor strengths of the two components) should occur in the 
case of a homogenous wind. The EW ratio of the Siiv 
doublet, which is observed to be smaller than 2, could be 
explained by the prese nce of optically thick structures 
present in the wind (cf. iPrinia fc Massall2010r) . 


iSurlan et al.l (|2013ll model the effect of macroclumps by 
means of 3D Monte Carlo simulations and obtain similar 
results. A consequence of accounting for macroclumps is 
the need to introduce a further parameter, Tmac, which 
specifies the separation between the clumps. The non- 
LTE nature of the line forrnation makes a simultaneous 
prediction of Tmac, D, and M practically impossible. We 
are therefore forced to adopt a value for one of these 
parameters. Since not much is known about the geome¬ 
try of the clumps, we avoid prespecifying Tmac- Instead, 
we choose to adopt D = 10 as a compromise between 
D ^ 2 and D ~ 100. Mo t ivated by hydrodynamic s tud- 
ies (|Feldmeier et al.l[l997l : IRunacres fc Owockil[2002D . we 
assume that the clumping initiates at r = 1.1 i?* and 
grows to its maximum contrast of D = 10 at r ^ 10 i?*. 
Our results depend only weakly on the depth-dependence 
of D(r). 

Fig.[T0] portrays the two Civ and Siiv resonance dou¬ 
blets as observed at phase (j) = 0.83 (blue line). The 
black dashed line depicts the synthetic composite spec¬ 
trum without the inclusion of macroclumping in the pri¬ 
mary. The red dotted line depicts the same composite 
spectrum, but with Lmac = 0.5 i?* adopted for the pri¬ 
mary. The strong effect of macroclumping on the res¬ 
onance lines is evident. The Ha line, as well as the 
photospheric features, are hardly affected by the macro¬ 
clumping formalism. The value Lmac = 0.5 R* provides 
a fair compromise for most wind lines, but the Si iv res¬ 
onance lines in fact require larger values of the order of 
Lmac ~ R*- The analysis therefore suggests L > 0.5 i?*. 
This value should not be given too much significance, as 
macroclumps are treated only as a rough approximation 
here. Nevertheless, we accounted for the major effects 
expected to rise from both optically thick and optically 
thin clumps, so mass-loss rates are unlikely to be very dif¬ 
ferent than those derived here. Future variability studies 
should help to further constrain the amount of inhomoge- 
niety in the primary’s wind. 



Fig. 10.— The observed Siiv and Civ resonance doublets (left 
and right panels, respectively) at phase ip = 0.83 (blue line) com¬ 
pared to the composite synthetic spectrum calculated with the pa¬ 
rameters in Table[T] but without macroclumping in the primary 
(black dashed line), and with macroclumping using Lmac = 0.5 iJ* 
(dotted red line). 


5.2. Macroclumping 

In PoWR, macroclumping is implemented only in the 
formal integration, leading to an effective reduct ion of 
the opacity in strong lines. lOskinova et al.l (120071) thor¬ 
oughly discuss the method, and illustrate the significant 
effect of macroclumping in the O supergiant Q Puppis. 


After assuming a clumping factor, the mass-loss rate 
is fairly well constrained, but still depends on the 
adopted distance. For the adopted distance of d = 
380 pc, the mass-loss rate is found to be logM = 
—6.4 ± 0.15 [Mq yr“^]. Since the mass-loss rate scales 
as we can consider the two extreme alterna¬ 

tives for the clumping factor 3 < L> < 50 to set lower 
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and upper bounds for the mass-loss rate of the pri¬ 
mary: —6.2 < log Ml < —6.8 [MQyr“^]. Interestingly, 
IVink et HI (|200Clll predict Myink = —6.48 [Mq yr“^] for a 
star with the parameters of the primary (for d = 380 pc) 
based on hydrodynamic calculations, which is in very 
good agreement with our results for the adopted clump¬ 
ing factor of D = 10 (cf. Tabled]). However, we note 
that IVink et aTI (120001) performed their calculations for 
homogenous winds, and it is not clear whether their pre¬ 
dictions would re main the same in th e case of significant 
inhomogeneities (jMuiires et al.l[2MTI) . 

6. WHERE DO THE X-RAYS IN 5 ORI A COME FROM? 

So far, our analysis focused on the “cool” stellar 
wind. However, the cool wind alone cannot account 
for the observed X-rays in S Ori A. It is commonly be¬ 
lieved that X-ray emission in single stars originates in 
the wind d ue to instability of th e line-driving mecha¬ 
nism le.g.. iFeldmeier et alJ 119971) or via acoustic driv - 
ing from subsurface convection (|Cantiello et al.l [2009D . 
In binary systems, an excess of hard X-ray flux may 
originate from w ind-wind collisions (|Williams et al.lFi990l : 
ICorcora^l2003l) . In this section, we explore the proper¬ 
ties of the “hot” X-ray producing component. 

6.1. Auger ionization 

The presence of strong X-ray radiation in OB type 
stars was hypothesized prior t o the first direct X- r ay ob - 
servations in massive stars. iCassinelli &: Olsoiil (|1979f) 
were the first to suggest that the detection of UV reso¬ 
nance lines of high-ionization ions such as N v and O vi, 
which are observed in many O type stars, may in dicate 
that Auger ionization (|Meitnerlll922t lAiigeHll923[) plays 
an important role in their formation. Auger ionization 
occurs when yery energetic photons remoye an electron of 
an inner shell (usually the K-shell), ultimately resulting 
in a double ionization. The Auger effect on UV spectra 
has b een frequently detect ed in studies of OB type stars 
(e.g., lOskinoya et al.l[2?)TTI ). 

The effect of X-rays in the wind of the primary 
Aal is eyident by merely inspecting its UV spectrum. 
The presence of the strong UV resonance doublets O vi 
AA1032,1038 and N v AA1239,1243 cannot be repro¬ 
duced otherwise. Indeed, an X-ray luminosity of ~ 
1.4- lO^^ergs”^ is reported for S Ori A (Paper I), which 
corresponds to ~ 10“® ®’^ times the total bolometric lumi¬ 
nosity of the systerrF^. The question remains, howeyer, 
as to the origin of the X-ray radiation observed. 

There are several arguments which suggest that emis¬ 
sion from wind-wind collisions do not dominate the X- 
ray flux in 6 Ori A. Paper H reports a relatively weak 
variability of the X-ray flux, which they cannot tie with 
certainty to phased-locked variations, in particular wind- 
wind collisions. The inferred value o f logEx/Unr,! 
—6.8 7 is typical for OB type s tars ( e. g.. iPallavicini et al l 
19811: iSeward fc Chlebowskil 119821 iMoffat et al.l 1200^ 
Oskinoval 120051 INaz^ I20m and does not imply any 
excess X-ray radiation from strong wind-wind collision. 
Lastly, the Nv AA1239,1243 resonance doublet, which 
clearly forms due to the X-ray radiation, is persistent 

We use here the total bolometric luminosity because all com¬ 
ponents are OB type stars and are thus expected to emit X-rays 
proportionally to their luminosities. 


throughout all available lUE spectra (see right panel of 
Fig. [TT] below). If most X-rays originated in a collision 
zone which irradiated only a part of the star (whose ori¬ 
entation relative to the observer depended on phase), 
then a larger variability could be anticipated. It there¬ 
fore seems plausible to assume that the X-rays originate 
in the wind itself. 

We model the eff ect of X-ray radiation as described by 
I Baum et al.l (I1992D . The X-ray emission is assumed to 
originate in optically thin filaments of shocked plasma 
embedded in the wind. The X-ray field is characterized 
by the temperature Tx of the shocked plasma and the 
radially-constant filling factor Agu, which describes the 
ratio of shocked to non-shocked plasma. The onset radius 
of shocked plasma is denoted by Rq. The X-ray emissiv- 
ity ?7x(»") at each radial layer r > Rq is proportional to 

. In principle, the three parameters Tx, Agp and Rq 
are chosen such that the observed X-ray SED (Huen- 
emoerder, priv. com.) is approximately reproduced by 
the synthetic X-ray SED emerging from the model, af¬ 
ter accounting for interstellar extinction. Once the onset 
radius Rq has been fixed, Tx and Agp follow from the 
observed SED. Motivated by the f/i analysis and X-ray 
line modeling results (see Sects. [631 and [S3|), we fix the 
onset radius to Rq = 1.1 i?*, leading to Tx = 3MK and 
Afiii = 0.1. These are only rough approximations to the 
X-ray properties in the wind (see Paper I), used to re¬ 
produce the bulk of X-ray emission observed. 



Fig. 11.— Left panel: Observed Nv resonance doublet at phase 
0 ~ 0.1 (blue line) with two s^thetic composite spectra, calculated 
with the parameters in Tabled not including X-rays (black dashed 
line), and including X-rays (red solid line). Right panel: Like the 
right panel of Fig.[7| the three lUE observations at phases (p = 
0.8,0.3, 0.1 (black, red, and green lines, respectively). 

The left panel of Fig.[TT] shows the effect of including 
X-rays in the wind of the primary on the N v resonance 
doublet AA1239,1243. The blue line depicts the lUE ob¬ 
servation of this doublet at phase (j) = 0.1. The black 
dashed line plots the synthetic composite spectrum with¬ 
out the inclusion of X-rays. Since the temperature in 
the wind is far from being sufficient to populate the N v 
ground state, only photospheric absorption is obtained. 
The red solid line is obtained after the inclusion of X- 
rays. While the line shape is not accurately reproduced, 
it is clear that X-rays are required in order to reproduce 
the observed P Cygni N v resonance line. Note also that 
the line is blended with components of iron lines in the 
blue part, making a determination of its terminal width 
difficult. The right panel of Fig.[TT] depicts three lUE 
observations at phases (j) ^ 0.8, 0.3, 0.1. While the Nv 
resonance doublet shows stronger variability than most 
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other resonance lines, the P Cygni feature is persistent 
and clearly visible in all 60 available lUE spectra. 

6 .2. f/i analysis 

Spectra of He-like ions (i.e. ions with two electrons) 
exhibit a group of three neighboring X-ray lines re¬ 
ferred to as resonance (r), forbidden (f), and intercom¬ 
bination (i) lines. A well established method to con¬ 
strain the formation region of X-rays in stellar winds 
makes use of observed forbidden-to-intercombination 
(f/i) line ratios, f/i analyses were or iginally developed 
to study the solar X-ray rad iation (|Gabriel fc Jordaiil 
119691: iBlumenthal et al.l Il972ll . but are now also fre¬ 
quently used to study t he X-ray emission of OB 
and Wolf-Rayet stars (e .g .. iWaldron fc Cassinelli [20011 : 
iLeutenegger et all 120061 IWaldron fc Cassinellil l20d7n . 
The f/i line ratio is determined by the relative popula¬ 
tion of the upper levels of the f and i lines, altered either 
by co llisions or by photo-excitations (e.g.. iPorauet et akl 
1200111 . For each helium-like ion, there are three possi¬ 
ble transitions between the upper levels of the f and i 
lines, denoted in the following with indices j = 0,1,2. 
These transitions are typically at UV frequencies, al¬ 
though their exact wavelengths depend on the ion. The 
stronger the UV radiation field is, the smaller the f/i 
ratio becomes. A detailed treatment allows one to con¬ 
struct an equation which predicts the f/i line ratio TZ as 
a function of t he radiative excitation r ate (j) and electron 
density rie (cf. IBlumenthal et al.lll972L Eq. Ic). 

The interpretation of observed f/i line ratios requires 
some discussion. After all, what we observe is likely 
the X-ray radiation reaching us from an extended region 
where it is formed. The simplest way to interpret the ob¬ 
servations would be to assume that the X-ray emitting 
gas is distributed over a thin spherical shell located at a 
formation radius inform, often referred to as the point-like 
interpretation. However, we know from detailed X-ray 
line fitting (see Sect 16.31) that the X-ray radiation must 
originate in an extended region. Obviously, the point-like 
interpretation cannot describe the whole truth. A more 
generalized interpretation to the observed f/ i line ra tios, 
thoroughly described by ILeutenegger et all (1200611 . in¬ 
volves the integration of the X-ray radiation emanating 
from a continuous range of radii. Within the assump¬ 
tions of this method, it is not the formation region which 
is sought, but the onset radius Rq o f the X-ray emis¬ 
sion. We refer the reader to stud ies byJGabriel_&_JordM 
( 196S 1. IBlumenthal et al.l (1197211 . and ILeutenegger et al.1 
(|2006ll for a detailed description of the methodology of 
both interpretations. 

We do not calculate the radiative excitation rate of 
the upper / level by dilu ting the photospheric fluxes, (cf. 
iLeuteiiegger et al1l2006L Eq. 2), but instead use directly 
the mean intensities at each radial layer, as obtained 
by our PoWR model. This way, we account for diffuse 
emission and limb darkening in a consistent manner. We 
include all three transitions j = 0,1, 2 in the calculation, 
properly weighted with their respective branching ratios. 
The relevant wavelengths Xj and oscillator strengths fj 
are extracted from the NIST database. 

Fig. [12] shows an example calculated for the He-like 
Mg XI ion. The blue dotted line depicts the f/i ratio 
TZ{r) predicted by the model for the point-like assump¬ 
tion as a function of the formation radius i?form, neglect- 



Fig. 12.— Theoretical f/i ratio 7^(r) is plotted as a function of 
the X-ray formation radius (blue dotted line) for the Mgxi ion. 
The intersection of 7^(r) with the observed value T^obs implies the 
formation radius under the assumptions of the point-like approach. 
The error Ur in the formation radius corresponds to l-cr measure¬ 
ment uncertainties in T^obs? depicted by the green shaded rectangle. 
The red solid line depicts 'JZ(r) as well, but includes the contribu¬ 
tion of collisions assuming a full ionization and a factor 1000 density 
enhancement in the shocked plasma 

ing collisional excitation. The solid horizontal line de- 
picts the observed value of TT-obs = 0.96 ± 0.36 (Paper 
I). The shaded green area depicts the l-cr measurement 
uncertainty. If the major part of the X-ray radiation 
originates at one radial layer, this layer will be located 
at i?form = where the uncertainty ar corre¬ 

sponds to the measurement uncertainty. In Fig. 1121 we 
also illustrate the influence of collisional excitation on 
TZ{r). The red solid line also plots TZ{r), but accounts 
for collisional excitation, assuming a full ionization and 
an unrealistically large factor of 1000 for the density en¬ 
hancement in the shocked regions. Evidently, the contri¬ 
bution of collisions to the excitation of the upper / level 
is negligible in S Ori A. 

In Eig. lI3l we compare the point-like interpretation 
with the more generaliz ed interpretation described by 
ILeutenegger et all (|2006D for the ion Mg xi. The red solid 
line depicts TZ{r), as in Eig.[T21 The green dashed line 
shows the predicted f/i ratio 7Z{r) as a function of the 
onset radius Rq. Note that while both functions are de¬ 
picted in the same coordinate frame, the meaning of r is 
different. As is evident from the figure, at the observed 
value of TZ = 0.96, the onset radius Rq predicted by in¬ 
tegrating and weighting the contribution of all layers at 
r > Rq is approximately 2.4 i?*. 

Finally, Fig. |T4| graphically summarizes our results for 
the He-like ions Si xiii. Mg xi, Ne ix, and O vii using both 
the point-like and generalized interpretations. The mea¬ 
sured f/i ratio of Sxv does not provide any constraints, 
so this ion is not included. Measured f/i values and their 
uncertainties are given in Paper 1. Formation radii in¬ 
ferred for each line are indicated with red diamonds, 
while onset radii are indicated with blue circles. The 
dashed vertical lines mark the 1 — a error range corre¬ 
sponding to f/i measurement uncertainties. The shaded 
gray area in both panels depicts the optically thick region 
in the model, i.e. regions where photons cannot escape. 

Evidently, the X-ray formation radii are different 
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region of X-ray formation (green dashed line). The intersection of 
'R-{r) with the observed value delivers the onset radius Rq. 

for ions with significantly different ionization potentials 
(e.g., Sixiii vs. Ovii). This may suggest that a distribu- 
tion of temperature s governs the forma tion radii (see e.g., 
iKrticka et al.ll20d^ iHerve et al.ll2013l ). Fig. [14] suggests 
that higher ions are formed at lower radii. This could 
imply that the hotter plasma is found closer to the stel¬ 
lar surface, while cooler plasma do minates farther out, a 
fact w hich was also pointed out bv I Waldron fc Cassinellil 
(|2007l) from their study of a large sample of 0-type X-ray 
spectra. The correlation between the ionization poten¬ 
tial and location in the wind is no longer seen for the 
onset radii, which may suggest that X-rays are emit¬ 
ted in all He-like lines already very close to the stellar 
surface, independent of the element. This is in agree¬ 
ment with a picture based on hydrodynamic simulations, 
where plasma of different temp eratures is present in c ool- 
ing layers behind a shock front (|Feldmeier et al.lll997fl . It 
is also worth noting that the formation radii roughly fol¬ 
low the same trend as the t = 1 surface of the cool wind, 
providing an independent confirmation that our analysis 
provides realistic values for the cool wind opacity. 

6.3. X-ray line modeling 

We now demonstrate that the high-resolution X-ray 
spectrum of 5 Ori A can be described consistently with 
the cool stellar wind model obtained from the analysis 
of UV and optical spectra. We use the high signal- 
to-noise ratio X-ray spectra obtained by the Chandra's 
HEG and MEG detectors. The observation consists of 
four different exposures taken at different orbital phases. 
A complete description of the X-ray observations and 
data analysis is given in Paper I. During the orbital mo¬ 
tion, the radial velocity of the primary changes between 
-|-114kms“^ and —78kms“^, with a systemic velocity of 
7 = 15kms“^ (Paper III). It is ambiguous whether the 
observed X-ray Doppler shifts follow the expected orbital 
pattern (Paper II). Regardless, our tests show that the 
resulting Doppler shifts, which are comparable to the 
instrumental resolution and are much smaller than the 
wind velocity, bear a negligible effect on the modeled line 
profiles. We therefore neglect Doppler shifts due to or¬ 
bital motion, and model the line profiles at the systemic 
velocity of 15kms“^. For all La lines of H-like ions, 



Fig. 14.— Inferred formation radii (red diamonds) and onset radii 
(blue circles) of the X-ray radiation, as inferred from the four He- 
like ions Si XIII, Mgxi, Neix, and O vii. The error bars correspond 
to measureme nt u ncertainties propagated into errors of the inferred 
radii (see Fig. ll3ll . The shaded area corresponds to radii at which 
"^Ross ^ 1* 


both components of the doublet are taken into account 
in the model. 

We analyze the X-ray spectra by simulating X- 
ray lines using our 2 -D stochastic stellar wind code 
(jOskinova et al.l l2004h and comparing them with ob¬ 
served lines. As we are interested in the line profiles, the 
maximum of the synthetic lines is normalized to their 
observed maximum. Our model is based on the assump¬ 
tion of a two component medium: the “cool wind” and 
the “hot wind”. The X-ray radiation is assumed to orig¬ 
inate in optically thin filaments evenly distributed in the 
hot wind component (the treat ment of resonant scatter - 
ing in X-ray lines is presented in lignace fc Gavle'^l2002t) . 
The emissivity ?7x scales as and the filling factor is 
assumed to be radially constant. A more sophisticated 
assumption on the behavior of px could lead to slightly 
different results numerically, but should not affect our 
conclusions qualitatively. The X-rays are attenuated on 
their way to the observer due to K-shell absorption in the 
cool wind. Hence the X-ray propagation in the stellar 
wind is described by a pure absorption case of ra diative 
transfer and is therefore r elatively simple (e.g., lignacel 
l200lHOskinova et al.ll200^ . 

iMacfarlane et ah ( 19911) demonstrated that the shape 
of optically thin emission lines is sensitive to the column 
density of the cool absorbing material: the lines become 
more skewed to the blue when the cool wind column den¬ 
sity is large. They suggested using this property to mea¬ 
sure the cool wind density. However, in clumped winds, 
the w ind column density is reduced (jEeldmeier et al.l 
Hool), implying that the analysis of observed X-ray emis¬ 
sion line profiles in 0-type stars should account for wind 
clumping. 

Several previous studies of X-ray emission lines in the 
spectra of 0-type stars reported that accounting for wind 
clumping in the models does not im prove the quality o f 
the fit to the X-ray spectrum (e.g., iHerve et al.l 120131) . 
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This is likely due to a degeneracy between the effects 
of mass-loss and clumping on the X-ray spectrum: one 
can always neglect clumping, albeit at the cost of lower 
estimated wind opacity. However, as we show here, the 
effect of clumping cannot be assumed to be negligible 
for the realistic mass-loss rate derived from the detailed 
analysis of the cool wind. 

In this study, we extract the mass-absorption coeffi¬ 
cient K\{r) directly from the PoWR model of the pri¬ 
mary’s cool wind. Importantly, our models include 
clumping and the feedback of X-ray radiation on the 
ionization structure of the wind and are therefore con¬ 
sistent. We make a simplifying assumption that the ve¬ 
locity field of the hot X-ray emitting plasma is the same 
as that of the cool X-ray absorbing wind. Thus, with the 
exception of the radial distribution of the hot plasma, all 
X-ray model parameters are determined by the non-LTE 
stellar atmosphere model. 

The theory of X-ray transfer in clumped s tellar 
winds was dev e loped by iFeldmeier et al.l (1200311 and 
lOskinova et al.l (1200411 . The special case of sph erical 
clumps was considered bv lOwocki fc Coh^ (|2006ll . The 
macroclumping formalism does not make any ad hoc 
assumptions about the size and optical depths of the 
cool wind clumps and is therefore suitable for effec¬ 
tively describing both smooth winds as well as clumped 
winds. This formalism further allows one to treat angle- 
dependent opacities, e.g., non-spherical wind clumps. A 
key parameter of the macroclumping formalism is the 
fragmentation frequency uq [s“^] - the number of clumps 
per unit time passing through a sphere at an arbitrary ra¬ 
dius, which does not depend on the radius due to the as¬ 
sumption of clump conservation, uq is approximately re¬ 
lated to the average separation between clumps Tmac (see 
Sect.[5|) via T^ac ~ Allowing for an angular- 

dependent opacity, e.g., flattened clumps or shell frag¬ 
ments, results in distinct line shapes, with flat topped 
line profiles in the limi ting case of opaque clumps (see 
also [ignace et al.l l20I2tl . For each line, we consider a 
smooth model an d a clumped model with anisotropic 
wind opacity (cf. lOskinova et al.l 1200411 . assuming flat - 
tened, “pancake-like” clumps (iLenine fc Moffat 1199911 . 
Comparing the observations with a grid of inhomoge- 
nious models, we find that Uq = 8.6 x 10 “®s“ provides 
a good compromise for several X-ray lines. This roughly 
corresponds to a wind flow time of tflow = R*/voo ~ 1.5 h 
(no ~ Amac ~ 1 i?*, Consistent with the 

lower bound Amac ^0-5 implied from the UV and opti¬ 
cal analysis. 

The cool wind is virtually transparent at wavelengths 
shorter than 10 A (see Fig.[T3]in Sect. 16.2D . In this case, 
the line shape can provide vital information regarding the 
velocity field and density distribution of the hot plasma. 
The left panel of Fig. [15] shows the observed Mgxii line 
compared with a smooth-wind model (blue dotted line) 
and a clumped-wind model (red solid line). The Mgxii 
line has the highest S/N ratio in the HEG wavelength 
range, and since the HEG has a superior spectral reso¬ 
lution (0.012 A) compared to the MEG (0.023 A), it is 
especially useful for studying the detailed shape of line 
profiles. Fig. ITS] illustrates the above statement: lines at 
A lOA are hardly absorbed in the wind, and thus do 
not show any sensitivity to wind clumping. Furthermore, 



Wavelength [A] Wavelength [A] 


Fig. 15.— Left panel: Observed co-added HEG±1 spectrum of 
(5Ori A centered on the Mgxii line (histogram), compared with 
a smooth-wind model (blue dotted line), and a clumped model 
(red solid line). As anticipated, clumping hardly affects the line 
formation at short wavelengths. Right panel: Again Mgxii, with 
three smooth models which assume X-ray onset radii of 1.01 U* 
(blue dotted-solid line), 1.1 fj* (red solid line) and 1.5 it* (green 
dashed line), suggesting that X-ray radiation initiates roughly at 
1.1 it* (for the adopted velocity law). The vertical lines indicate 
the rest wavelength shifted with the systemic velocity 15kms“^ 

the good agreement between the models and the obser¬ 
vations implies that the majority of X-rays originate in 
filaments co-moving with the cold wind, exhibiting the 
whole range of velocities up to the roughly the terminal 
velocity in their line profiles. 

The right panel of Fig.[T5|shows a comparison between 
the observed Mgxii line and three smooth-wind models 
that assume different onset radii for the X-ray emission: 
1.01 A* (blue dotted solid line), 1.1 i?* (red solid line), 
and 1.5 i?* (green dashed line). This comparison clearly 
illustrates that X-rays must form already very close to 
the stellar surface, at around 1.1 i?*. Larger onset radii 
result in flat-topped line profiles which do not agree at 
all with the observed profiles. Smaller onset radii imply 
narrower lines than those observed. We note that the 
exact values are very sensitive to the adopted velocity 
law. However, the same conclusion is obtained for all 
our test models, where we vary the exponent of the /?- 
law in the domain 0.5 < /3 < 1.5: X-rays are formed in 
the wind, and onset radii are close to the stellar surface. 
Note that a small onset radius does not imply that the 
X-ray forms only close to the stellar surface. 



Wavelength [A] Wavelength [A] 


Fig. 16.— Left panel: Observed co-added MEG±1 spectrum of 
5 Ori A centered on Fe xvil line (histogram) and two wind model 
lines: smooth (blue dotted line) and clumped (red solid line). The 
onset of X-ray emitting plasma is at 1.2 it*. Right panel: Same as 
left panel, but for the O VIII line 

While the effect of porosity is negligible at A < 
10 A, it should generally be taken into account when 
modeling lines at longer wavelengths, where the cool 
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wind opacity increases. Fig.[TH] shows a comparison of 
the observed Fexvii A15.01A (left panel) and Oviii 
AA18.967,18.972A ( right panel) lines with models that 
include and neglect clumping (red solid and blue dot¬ 
ted lines, respectively). The wind model that neglects 
clumping is in poor agreement with the data, while a 
better agreement is reached when clumping is assumed. 
Nevertheless, the exact profile shape is not reproduced, 
especially in the case of the O viii line. In fact, there 
is evidence that the widths of the X-ray lines in S Ori A 
change with time (Paper II). Our tests indicate that the 
observed lines can be more accurately reproduced when 
htted at each phase separately. To fit the global pro- 
hles, the variable line profiles would need to be averaged 
out. However, such a study, while interesting, is not the 
subject of this Paper. 

To summarize, we illustrated that clumping generally 
plays a role in shaping the observed X-ray lines, that the 
line profiles suggest small onset radii of ^ 1.1 i?*, and 
that the X-ray emitting parcels are likely coupled to the 
cool wind. 

7. SUMMARY AND CONCLUSIONS 

We have performed a multi-wavelength, non-LTE spec¬ 
troscopic analysis of the massive d Ori A system, which 
contains the visually brightest 0-type eclipsing binary in 
the sky and a wider tertiary component. Our goal was 
to obtain accurate stellar parameters for the components 
of the system, to analyze their winds, and to gain a bet¬ 
ter understanding of the X-ray radiation emitted from 
the system. Three additional studies performed within 
the framework of the current collaboration explore the 
X-ray properties (Paper I) and variability (Paper II) of 
the system, and conduct a complete binary and optical 
variabilty analysis of the system (Paper III). 

We conclude the following: 

• Distance-independent parameters such as T* and 
log g^f derived for the primary (cf. Table[T]) are con¬ 
sistent with its spectral type 09.5 II. Distance- 
dependent parameters such as the luminosity and 
mass are found to be unusually low for the Hip- 
parcos distance d = 212 pc. Typical values (Li = 
5.28 Lq, Ml = 24 Mq) are obtained for the dis¬ 
tance d = 380 pc of the neighboring cluster. These 
results agree well with the results reported in Pa¬ 
per III, suggesting that the Hipparcos distance is 
strongly underestimated. 

• The secondary is marginally detectable in the com¬ 
posite spectrum. The V-band magnitude differ¬ 
ence between the primary and secondary is con¬ 
strained to AVAaiAa 2 2™8. The parameters of 
the secondary suggest it is a B type dwarf. The 
tertiary is confirmed to be a rapid rotator with 
usini ~ 220kms“^, and its parameters corre¬ 
spond to an early B type or a late O type subgiant. 

• Rapid turbulent velocities (~ 200kms“^) prevail 
in the wind of the primary close to the stellar sur¬ 
face. We further find evidence for optically thick 
wind-inhomogeneities (“macroclumping”), affect¬ 
ing both strong resonance lines and X-ray lines. 


• For a clumping factor of D = 10 and accounting 
for porosity, the primary’s inferred mass-loss rate 
is logM ft! —6.4 [Mq yr“^]. This value provides 
a consistent picture for S Ori A along all spectral 
domains, from X-rays and UV to the optical and 
radio, and is furthermore in good agreement with 
hydrodynamical predictions. 

• Most X-rays emerging from S Ori A are likely in¬ 
trinsic to the wind of the primary. X-ray onset radii 
are found to be ^ l.li?*. 

While 6 Ori A gives us a precious opportunity to study 
stellar winds and multiple systems of massive stars, there 
are still questions left unanswered. Does the primary 
really have an exceptionally low mass, luminosity, and 
mass-loss rate, or is the Hipparcos distance significantly 
underestimated? Did mass-transfer occur between the 
primary Aal and the secondary Aa2? Where do the 
clumps in the primary wind originate, and how can they 
be further quantified? Why do we not clearly observe X- 
ray emission originating in the interaction between the 
primary and secondary? How did the system evolve to 
its current state, and how will it continue to evolve? To 
answer these questions, efforts should be made in obtain¬ 
ing high S/N optical and X-ray spectra which resolve 
the binary Aa from the tertiary Ab, and which enable a 
proper disentanglement of the system. Moreover, further 
studies should provide more hints regarding the correct 
distance of the system. Finally, variability studies of P 
Cygni lines in the system should provide crucial, inde¬ 
pendent constraints on the amount of inhomogeniety in 
the primary’s wind. Answering these questions should 
take us one step closer to a more complete understand¬ 
ing of stellar winds, and of the evolution and properties 
of massive stars in multiple systems. 
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